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Quantum Error Correction and fault-tolerant quantum computation represent arguably the most 
vital theoretical aspect of quantum information processing. It was well known from the early 
developments of this exciting field that the fragility of coherent quantum systems would be a 
catastrophic obstacle to the development of large scale quantum computers. The introduction of 
Quantum Error Correction in 1995 showed that active techniques could be employed to mitigate 
this fatal problem. However, Quantum Error Correction and fault-tolerant computation is now a 
much more mature field and many new codes, techniques and methodologies have been developed 
to implement error correction for large scale quantum algorithms. This development has been so 
pronounced that many in the field of quantum information, specifically those new to quantum 
information or those focused on the many other important issues in quantum computation have 
not been able to keep up with the general formalisms and methodologies employed in this area. 
In response we have attempted to summarize the basic aspects of Quantum Error Correction 
and fault-tolerance, not as a detailed guide, but rather as a basic introduction. Rather than 
introducing these concepts from a rigorous mathematical and computer science framework, we 
instead examine error correction and fault-tolerance largely through detailed examples, progressing 
from basic examples of the 3-qubit code through to the stabilizer formalism which is now extremely 
important when understanding large scale concatenated code structures, quantum circuit synthesis 
and the more recent developments in subsystem and topological codes. 
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The micro-computer revolution of the late 20th century 
has arguably been of greater impact to the world that any 
other technological revolution in history. The advent of 
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transistors, integrated circuits, and the modern micro- 
processor has spawned literally hundreds of devices from 
pocket calculators to the iPod, all now integrated through 
an extensive worldwide communications system. How- 
ever, as we enter the 21st century, the rate at which com- 
putational power is increasing is driving us very quickly 
to the realm of quantum physics. The component size 
of individual transistors on modern microprocessors are 
becoming so small that quantum effects will soon begin 
to dominate over classical electronic properties. Unfor- 
tunately the current designs for micro-electronics mean 
that quantum mechanical behavior will tend to result 
in unpredictable and unwanted behavior. Therefore, we 
have two choices: to keep trying to suppressing quan- 
tum effects in classically fabricated electronics or move 
to the field of quantum information processing (QIP) 
where we instead exploit them. This leads to a paradigm 
shift in the way we view and process information and 
has lead to considerable interest from physicists, en- 
gineers, computer scientists and mathematicians. The 
counter-intuitive and strange rules of quantum physics 
offers enormous possibilities for information processing 
and the development of a large scale quantum computer 
is the holy grail of many groups worldwide. 

While the advent of Shor's algorithm (Sho97) certainly 
spawned great interest in quantum information process- 
ing and demonstrated that the utilization of a quan- 
tum computer could lead to algorithms far more efficient 
than those used in classical computing, there was a great 
deal of debate surrounding the practicality of building a 
large scale, controllable, quantum system. It was well 
known even before the introduction of quantum informa- 
tion that coherent quantum states were extremely fragile 
and many believed that to maintain large, multi-qubit, 
coherent quantum states for a long enough time to com- 
plete any quantum algorithm was unrealistic (Unr95). 
Additionally, classical error correction techniques are in- 
trinsically based on a digital framework. Hence, can the 
vast amount of knowledge gained from classical coding 
theory be adapted to the quantum regime where while 
the readout of qubits is digital but actual manipulations 
are analogue. 

Starting in 1995, several papers appeared, in rapid suc- 
cession, proposing codes which were appropriate to per- 
form error correction on quantum data (Sho95; Ste96a; 
ICS96I ILMPZ96P - This was the last theoretical aspect 
needed to convince the general community that quantum 
computation was indeed a possibility. Since this initial 
introduction, the progress in this field has been extensive. 

Initial work on error correction focused heavily on de- 
veloping quantum codes (|Ste96bl ICG97I IGot96l IPVK97[) , 
introducing a more rigorous theoretical framework for 
the structure and properties of Quantum Error Correc- 
tion (QEC) (iKLOOl ICRSS98I IGro97l IKLV001 IKLP05I) 
and the introduction of concepts such as fault-tolerant 
quantum computation (|Sho96| IDS96| IGot98p which 
leads directly to the threshold theorem for concatenated 
QEC (|KLZ96irATI097| . In more recent years QEC pro- 



tocols have been developed for various s ystems, su ch as 
continuous variables (|LS981 IBraMl lvL08l |ATK+08|) . ion- 
traps and other systems containing motional degrees of 
freedom (|LW03bl |ST98| . adiabatic computatio n (jJFS06[) 
and globally controlled quantum computers (|BBK03[1 . 
Additionally, work still continues on not only developing 
more complicated (and in some ways, more technologi- 
cally useful) protocols such as subsystem codes ()Bac06[) 
and topological codes (IKit97l IDKLP02I IRHG071 IFSG08I) 
but also advanced techniques to implement error correc- 
tion in a fault-tolerant manner (|Ste97a| ISte02> IDA07jl . 

Along with QEC, other methods of protecting quan- 
tum information were also developed. These other tech- 
niques would technically be placed in a separate cate- 
gory of error avoidance rather than error correction. The 
most well known technique of error avoidance is proto- 
cols such as decoherence free subspaces (DFS) (DG97 
IDG98bl IZR97bl IZR97al IDG98al IIW03a[) . While this 
protocol has the mathematical structure of a self correct- 
ing quantum code, it is largely a technique to suppress 
certain, well structured, noise models. As with QEC, this 
field of error avoidance is vast, now incorporating ideas 
from optimal control to create specially designed control 
sequences to counteract the effect of errors induced from 
environmental coupling. These new methods of dynami- 
cal decoupling can take simple structures such as Bang- 
Bang control (|VL98I IVT99IIZan99|l . to more complicated 
and generalized protocols to help decouple qubits from 
the environment (|VKL99I IFLP04I IVK03I IVK05|) . 

This review deals exclusively with the concepts of 
QEC and fault-tolerant quantum computation. Many 
papers have reviewed error corre ction and Fault- 
tolerance ()Got97l INC00I IGotf)2l lKLA+021 ISteOTl 
IGot09p . however to cater for a large audience, we at- 
tempt to describe QEC and Fault-tolerance in a much 
more basic manner, largely through examples. Instead 
of providing a more rigorous review of error correction, 
we instead try to focus on more practical issues involved 
when working with these ideas. For those who have re- 
cently begun investigating quantum information process- 
ing or those who are focused on other important theo- 
retical and/or experimental aspects related to quantum 
computing, searching through this enormous collection of 
work is daunting especially if a basic working knowledge 
of QEC is all that is required. We hope that this review 
of the basic aspects of QEC and fault-tolerance will allow 
those with little knowledge of the field to quickly become 
accustomed to the various techniques and tricks that are 
commonly used. 

We begin the discussion in section [TT] where we share 
some preliminary thoughts on the required properties of 
any quantum error correcting protocol. In section III 
we review some basic noise models from the context of 
how they influence quantum algorithms. Section |VII| 
introduces quantum error correction through the tradi- 
tional example of the 3-qubit code, illustrating the cir- 
cuits used for encoding and correction and why the prin- 
cipal of redundant encoding suppresses the failure of en- 
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coded qubits. Section |"VII| then introduces the stabilizer 
formalism (Got97), demonstrating how QEC circuits are 
synthesized once the structure of the code is known. In 
section |IX| we then briefly return to the noise models 
and relate the abstract analysis of QEC, where errors 
are assumed to be discrete and probabilistic, to some of 
the physical mechanisms which can cause errors. Sec- 
tions [X] and |XI| introduces the concept of fault-tolerant 
error correction, the threshold theorem and how logi- 
cal gate operations can be applied directly to quantum 
data. We then move on to circuit synthesis in section |XTT| 
presenting a basic fault-tolerant circuit design for log- 
ical state preparation using the [[7,1,3]] Steane code 
as a representative example of how to synthesize fault- 
tolerant circuits from the stabilizer structure of quan- 
tum codes. Finally in section |XIV| we review specific 
codes for qubit loss and examine two of the more mod- 
ern techniques for error correction. We briefly examine 
quantum subsystem codes (|Bac06p and topological sur- 
face codes (|DKLP02t IFSG08[) due to both their theoret- 
ical elegance and th eir increas ing relevance in quantum 
architecture designs (|DFS + 08|l . 



II. PRELIMINARIES 

Before discussing specifically the effect of errors and 
the basics of Quantum Error Correction (QEC) we first 
dispense with the very basics of qubits and quantum 
gates. We assume a basic working knowledge with quan- 
tum information (EJ96 ; NCOOi) and this brief discussion 
is used simply to define our notation for the remainder 
of this review. 

The fundamental unit of quantum information, the 
qubit, which unlike classical bits can exist in coherent su- 
perpositions of two states, denoted |0) and |1). These ba- 
sis states can be photonic polarization, spin states, elec- 
tronic states of an ion or charge states of superconducting 
systems. An arbitrary state of an individual qubit, |</>), 
can be expressed as, 

\4,} = a\O)+0\l) (1) 

where normalization requires, |a| 2 + |/3| 2 — 1. Quan- 
tum gate operations are represented by unitary opera- 
tions acting on the Hilbert space of a qubit array. Unlike 
classical information processing, conservation of proba- 
bility for quantum states require that all operations be 
reversible and hence unitary. 

When describing a quantum gate on an individual 
qubit, any dynamical operation, G, is a member of the 
unitary group U(2), which consists of all 2 x 2 matrices 
where G^ = G~ l . Up to a global (and unphysical) phase 
factor, any single qubit operation can be expressed as a 
linear combination of the generators of SU(2) as, 

G = C/(T/ + C X (T X + CylJy + C Z (T Z (2) 



where, 

°* = (l J) ' a v = \i o) ' ff » = (J -l) ' (3) 

are the Pauli matrices, 07 is the 2x2 identity matrix and 
the co-efficients (ci,c x ,c y ,c z ) £ C satisfy |c/| 2 + Ic^l 2 + 

\Cy\ 2 + \ Cz \ 2 = l. 

The concept of Quantum Error Correction (QEC) is 
fundamental to the large scale viability of quantum in- 
formation processing. Although the field is largely based 
on classical coding theory, there are several issues that 
need to be considered when transferring classical error 
correction to the quantum regime. 

First, coding based on data-copying, which is exten- 
sively used in classical error correction cannot be used 
due to the no-cloning theorem of quantum mechan- 
ics (jWZ82 ) . This result implies that there exists no trans- 
formation resulting in the following mapping, 

0»®M = |0)®|0>. ( 4 ) 

i.e. it is impossible to perfectly copy an unknown quan- 
tum state. This means that quantum data cannot be 
protected from errors by simply making multiple copies. 
Secondly, direct measurement cannot be used to effec- 
tively protect against errors, since this will act to de- 
stroy any quantum superposition that is being used for 
computation. Error correction protocols must therefore 
be employed which can detect and correct errors without 
determining any information regarding the qubit state. 
Finally, unlike classical information, qubits can experi- 
ence traditional bit errors, |0) <-> |1) but unlike classical 
information, qubit are also susceptible to phase errors 
|1) <-> — Hence any error correcting procedure needs 
to be able to simultaneously correct for both. 

At its most basic level, QEC utilizes the idea of redun- 
dant encoding where quantum data is protected by ex- 
tending the size of the Hilbert space for a single, logically 
encoded qubit and essentially spreading out the informa- 
tion over multiple qubits. This way, errors only perturb 
codeword states by small amounts which can then be 
detected and corrected, without directly measuring the 
quantum state of any qubit. 

III. QUANTUM ERRORS: CAUSE AND EFFECT 

Before we even begin discussing the details of quantum 
error correction, we first examine some of the common 
sources of errors in quantum information processing and 
contextualize what they imply for computation. We will 
consider several important sources of errors and how they 
influence a trivial, single qubit, quantum algorithm. 

This trivial algorithm will be a computation consisting 
of a single qubit, intitilized in the |0) state undergoing N 
identity operations. Such that the final, error free state 
is, 

N 

i^) fina i = i°) = i°> ' ( g ) 
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Measurement of the qubit in the |0), |1) basis will conse- 
quently yield the result with a probability of unity. We 
examine, independently, several common sources of error 
from the effect they have on this simple quantum algo- 
rithm. Hopefully, this introductory section will show that 
while quantum errors are complicated physical effects, in 
QIP the relevant measure is the theoretical success prob- 
ability of a given quantum algorithm. 

A. Coherent Quantum Errors: You don't know what you 
are doing! 

The first possible source of error is coherent, systematic 
control errors. This type of error is typically associated 
with bad system control and/or characterization where 
imprecise manipulation of the qubit introduces inaccu- 
rate Hamiltonian dynamics. As this source of error is 
produced by inaccurate control of the system dynamics 
it does not produce mixed states from pure states (i.e. 
it is a coherent, unitary error and does not destroy the 
quantum coherence of the qubit but instead causes you 
to apply an undesired gate operation). In our trivial al- 
gorithm, we are able to model this several different ways. 
To keep things simple, we assume that incorrect char- 
acterization of the control dynamics leads to an identity 
gate which is not 07, but instead introduces a small ro- 
tation around the X-axis of the Bloch sphere, i.e. 

N 

Mtoai = IP"* 1°) = °°< Ne ) 1°) + ^in(iVe) |l) . (6) 

We now measure the system in the |0) or |1) state. In 
the ideal case, the computer should collapse to the state 
|0) with a probability of one, P(|0)) = 1. However we 
now find, 

P(|0)) = cos 2 (7V e )«l~(iV e ) 2 , 

F(|l)) = sin 2 (iVe) « (iVe) 2 . 

Hence, the probability of error in this trivial quantum 
algorithm is given by p e rror ~ (Ne) 2 , which will be small 
given that Ne -C 1. The systematic error in this system 
is proportional to both the small systematic over rotation 
and the total number of applied identity operations. 

B. Decoherence: The devil is in the environment 

Environmental decoherence is another important 
source of errors in quantum systems. Once again we will 
take a very basic example of a decoherence model and 
examine how it influences our trivial algorithm. Later in 
section |IX| we will illustrate a more complicated decoher- 
ence model that arises from standard mechanisms. 

Consider a very simple environment, which is another 
two level quantum system. This environment has two 
basis states, |eo) and \e\) which satisfies the completeness 
relations, 

(e i \e j ) = 8 ij , |e )(e | + |e 1 )(e 1 |=/. (8) 



We will also assume that the environment couples to the 
qubit in a specific way. When the qubit is in the |1) 
state, the coupling flips the environmental state while if 
the qubit is in the |0) state nothing happens to the envi- 
ronment. Additionally, as we anticipate the effect of this 
decoherence model we will slightly alter our trivial algo- 
rithm. Rather than considering a qubit prepared in the 
|0) state and applying N identity operations, we instead 
modify the algorithm to the following, 

IV>> fina i = HajH |0) = Haj^-dO) + |1)) 

(9) 

= #-L(|0) + |l» = |0>. 

Essentially we are performing two H = Hadamard op- 
erations separated by a wait stage, represented by the 
Identity gate. Finally, this model assumes the sys- 
tem/environment interaction only occurs during this wait 
stage of the algorithm. As with the previous algorithm we 
should measure the state |0) with probability one. The 
reason for modifying our trivial algorithm is because this 
specific decoherence model acts to reduce coherence be- 
tween the |0) and |1) basis states and hence we require 
a coherent superposition to observe any effect from the 
environmental coupling. 

We now assume that the environment starts in the pure 
state, \E) — |eo), and couples to the system such that, 

Ha r H |0) \E) = i(|0) + |1» |e ) + ~(|0) - |1» | ei ) (10) 

As we are considering environmental decoherence, pure 
states will be transformed into classical mixtures, hence 
we now move into the density matrix representation for 
the state HaiH\0) \E), 

Pf=~W) (0| + |0)(1| + |1) (0| + |l){l|)|e ) (eo| 

+ i(|0) (0|-|0)(l|-|l)(0| + |l)(l|)| ei )( ei | 

f (11) 
+ i(|0) (0|-|0)(l| + |l)(Q|-|l)(l|)|e ) ( ei | 

+ i(|0) (0| + |0)(l|-|l)(0|-|l)(l|)|e 1 )(e o |. 

Since we do not measure the environmental degrees of 
freedom, we trace over this part of the system, giving, 

Tr E (p f ) = \(\0) (0| + |0>(1| + |1>(0| + |1)(1|) 

+ |(|0)(0|-|0)(1|-|1)(0| + |1)(1|) (12) 

= ^(|0) (0| + |1)(1|). 

Measurement of the system will consequently return |0) 
50% of the time and |1) 50% of the time. This final state 
is a complete mixture of the qubit states and is conse- 
quently a classical system. The coupling to the environ- 
ment removed all the coherence between the |0) and |1) 
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states and consequently the second Hadamard transform, 
intended to rotate (|0) + |l))/\/2 — > |0) has no effect. 

Since we assumed that the system/enviroment cou- 
pling during the wait stage causes the environmental de- 
gree of freedom to "flip" when the qubit is in the |1) 
state, this decoherence model implicitly incorporates a 
temporal effect. The temporal interval of our identity 
gate in the above algorithm is long enough to enact this 
full controlled- flip operation. If we assumed a controlled 
rotation that is not a full flip on the environment, the 
final mixture will not be 50/50. Instead there would 
be a residual coherence between the qubit states and an 
increased probability of our algorithm returning a |0). 
Section If XI revisits the decoherence model and illustrates 
how time-dependence is explicitly incorporated. 

C. Loss, Leakage, Measurement and Initialization: 
Variations of the above 

Other sources of error such as qubit initialization, mea- 
surement errors, qubit loss and qubit leakage are modeled 
in a very similar manner. Measurement errors and qubit 
loss are modeled in the same way as environmental deco- 
herence. Measurement errors are described by utilizing 
the following measurement projection onto a qubit space, 

A=(1- Pm )|0)(0|+pm|1)<1| (13) 

where pu € [0, 1] is the probability of measurement error. 
If we have a pure state p = |0) (0|, the probability of 
measuring a |0) is, 

P(\0))=Tr(Ap) = (l-p) (14) 

indicating that the correct result is observed with prob- 
ability 1 — p. 

Qubit loss is modeled in a slightly similar manner. 
When a qubit is lost, it is essentially coupled to the en- 
vironment which acts to measure the system, with the 
classical information lost. This coupling follows the de- 
coherence analysis shown earlier, where a 50/50 mixed 
state of the qubit results. Therefore the projector onto 
the qubit space is given by A — ^(\0) (0| + |1) (1|), which is 
identical to simply tracing over the lost qubit and equiv- 
alent to a measurement error of probability p = 0.5. 
With this type of error channel, not only is the physi- 
cal object lost (and hence cannot be directly measured), 
but an initially pure qubit is converted to a completely 
mixed state. While this model of qubit loss is equiva- 
lent to environmental coupling, correcting this type of 
error requires additional machinery on top of standard 
QEC protocols. The difficulty with qubit loss is the ini- 
tial detection of whether the qubit is actually present. 
While standard correction protocols can protect against 
the loss of information on a qubit, this still assumes that 
the physical object still exists in the computer. Hence 
in loss correction protocols, an initial non-demolition de- 
tection method must be employed (which determines if 



the qubit is actually present without performing a pro- 
jective measurement on the computational state) before 
standard correction can be utilized to correct the error. 

Initialization of the qubit can be modeled either using a 
coherent systematic error model or using the decoherence 
model. The specific methodology depends largely on the 
physical mechanisms used to initialize the system. If a 
decoherence model is employed, initialization is modeled 
exactly the same way as imperfect measurement. If we 
have a probability pi of initialization error, the initial 
state of the system is given by the mixture, 

p ! = (l- W )|0)(0|+p 7 |l)(l|. (15) 

In contrast, we could consider an initialization model 
which is achieved via a coherent unitary operation where 
the target is the desired initial state. In this case, the 
initial state is pure, but contains a non-zero amplitude 
of the undesired target, for example, 

|V^ = a|0)+/?|l) (16) 

where \&\ 2 + 1/?| 2 = 1 and \(3\ 2 <C 1. The interpretation of 
these two types of initialization models is identical to the 
coherent and incoherent models presented. Again, the 
effect of these types of errors relates to the probabilities 
of measuring the system in an erred state. 

One final type of error that we can briefly mention is 
the problem of qubit leakage. Qubit leakage manifests 
itself due to the fact that most systems utilized for qubit 
applications are not simple two level quantum systems. 
For example, Fig[l](from Ref. (|Ste97b[0 illustrates the en- 
ergy level structure for a 43 Ca + ion utilized for ion trap 
quantum computing at Oxford. The qubit in this sys- 
tem is defined only with two electronic states, however 
the system itself contains many more levels (including 
some which are used for qubit readout and initialization 
through optical pumping and photo- luminescence) . As 
with systematic errors, leakage can occur when improper 
control is applied to such a system. In the case of ion- 
traps, qubit transitions are performed by focusing finely 
tuned lasers resonant on the relevant transitions. If the 
laser frequency fluctuates or additional levels are not suf- 
ficiently detuned from the qubit resonance, the following 
transformation could occur, 

CT|0) = a|0)+/3|l)+ 7 |2), (17) 

where the state |2) is a third level which is now populated 
due to improper control. The actual effect of this type of 
error can manifest in several different ways. The primary 
problem with leakage is that it violates the basic assump- 
tion of a qubit structure to the computer. As quantum 
circuits and algorithms are fundamentally designed as- 
suming the computational array is a collection of 2-level 
systems, operators of the above form (which in this case 
is operating over a 3-level space) will naturally induce 
unwanted dynamics. Another important implication of 
applying non-qubit operations is how these levels inter- 
act with the environment and hence how decoherence 
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1 (from Ref. (|Ste97b[> ) Energy level structure for the 
Ca + investigated by the Oxford ion-trapping group. The 
structure of this ion is clearly not a 2-level quantum system. 
Hence leakage into non-qubit states is an important factor to 
consider. 



effects the system. For example, in the above case, the 
unwanted level, |2), may be extremely short lived leading 
to an emission of a photon and the system relaxing back 
to the ground state. For these reasons, leakage is one 
of the most problematic error channels to correct using 
QEC. In general, leakage induced errors need to be cor- 
rected via the non-demolition detection of a leakage event 
(i.e. determining if the quantum system is confined to 
a qubit without performing a measurement discriminat- 
ing the |0) and |1) states (|Pre98l IGBP97I IVWW05|) ') or 
through the use of complicated pulse control which acts 
to re-focus a improperly confined quantum gate back to 
the qubit subspace (IWBL02|lBLWZ05p . In the context of 
mass manufacturing of qubit systems, leakage would be 
quantified immediately after the fabrication of a device, 
using intrinsic cha racterizat ion protocols such as those 
discussed in Ref. (|DSO + 07|) . If a particular system is 
found to be improperly confined to the qubit subspace it 



would simply be discarded. Employing characterization 
at this stage would then eliminate the need to implement 
pulse control of leakage, shortening gate times and ulti- 
mately reducing error rates in the computer. 

In this section we introduced the basic ideas of quan- 
tum errors and how they effect the success of a quantum 
algorithm. Section |LX| will return in a more focused man- 
ner to error models and how they relate to error correc- 
tion in a quantum computer. 



IV. QEC, 
CODE 



A GOOD STARTING POINT: THE 3-QUBIT 



The 3-qubit bit-flip code is traditionally used as a basic 
introduction to the concept of Quantum Error Correc- 
tion. However, it should be emphasized that the 3-qubit 
code does not represent a full quantum code. This is due 
to the fact that the code cannot simultaneously correct 
for both bit and phase flips (see section. IX), which is a 
sufficient condition for correcting errors for an arbitrary 
error mapping on a single qubit. This code is a standard 
repetition code which was extended by Shor (ISho95jl to 
the full 9-qubit quantum code which was the first demon- 
stration that QEC was possible. 

The 3-qubit code encodes a single logical qubit into 
three physical qubits with the property that it can correct 
for a single a x = X bit-flip error. The two logical basis 
states |0) L and |1) L are defined as, 



= 1000) , 



11^ = 1111) 



such that an arbitrary single qubit state \% 
is mapped to, 

a\0) + (3\l)^a\0) L + (3\l) L 

= a|000)+/3|lll) = 



(18) 
a\Q)+0\l) 



(19) 



Fig. [2] illustrates the quantum circuit required to encode 
a single logical qubit via the initialization of two an- 
cilla qubits and two CNOT gates. The reason why this 
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FIG. 2 Quantum Circuit to prepare the \Q) L state for the 3- 
qubit code where an arbitrary single qubit state, is coupled 
to two freshly initialized ancilla qubits via CNOT gates to 
prepare | ip) , . 



code is able to correct for a single bit flip error is the 
binary distance between the two codeword states. No- 
tice that three individual bit flips are required to take 



|0), 



1) L , hence if we assume \ip) = |0) L , a single bit 
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flip on any qubit leaves the final state closer to |0) L than 
|1) L . The distance between two codeword states, d, de- 
fines the number of errors that can be corrected, t, as, 
t — [(d — 1)/2J . In this case, d = 3, hence t = 1. 

How are we able to correct errors using this code with- 
out directly measuring or obtaining information about 
the logical state? Two additional ancilla qubits are in- 
troduced, which are used to extract syndrome informa- 
tion (information regarding possible errors) from the data 
block without discriminating the exact state of any qubit, 
Fig. [3] illustrates. For the sake of simplicity we assume 
that all gate operations are perfect and the only place 
where the qubits are susceptible to error is the region 
between encoding and correction. We will return to this 
issue in section |X] when we discuss Fault-tolerance. We 
also assume that at most, a single, complete bit flip er- 
ror occurs on one of the three data qubits. Correction 
proceeds by introducing two ancilla qubits and perform- 
ing a sequence of CNOT gates, which checks the parity 
of the three qubits. Table |T] summarizes the state of the 
whole system, for each possible error, just prior to mea- 
surement. 



Error Location 


Final State, |data) | ancilla) 


No Error 
Qubit 1 
Qubit 2 
Qubit 3 


a 000) 
a 100) 
a 010) 
a 001) 


00) +/3 
11) +0 
W) +(3 

01) + 


111) |00) 
011)|11) 
101) |10) 
110) 01) 



TABLE I Final state of the five qubit system prior to the 
syndrome measurement for no error or a single X error on 
one of the qubits. The last two qubits represent the state 
of the ancilla. Note that each possible error will result in a 
unique measurement result (syndrome) of the ancilla qubits. 
This allows for a X correction gate to be applied to the data 
block which is classically controlled from the syndrome result. 
At no point during correction do we learn anything about a 
or (3. 



For each possible situation, either no error or a single 
bit-flip error, the ancilla qubits are flipped to a unique 
state based on the parity of the data block. These qubits 
are then measured to obtain the classical syndrome re- 
sult. The result of the measurement will then dictate if 
an X correction gate needs to be applied to a specific 
qubit, i.e. 



Ancilla Measurement: 1 00) , 

Ancilla Measurement: |01) , 

Ancilla Measurement: 1 10) , 

Ancilla Measurement: 1 11) , 



Collapsed State 
Collapsed State 
Collapsed State 
Collapsed State 



a|000) +£1111) 
a|001) + /3|110) 
a|010) + /3|101) 
a 1 100) +/3|011) 



Clean State 
Bit Flip on Qubit 3 
Bit Flip on Qubit 2 
Bit Flip on Qubit 1 



(20) 



Provided that only a single error has occurred, the 
data block is restored. Notice that at no point dur- 
ing correction do we gain any information regarding the 
co-efficients a and f3, hence the computational wave- 
function will remain intact during correction. 

This code will only work if a maximum of one error 
occurs. If two X errors occur, then by tracking the circuit 
through you will see that the syndrome result becomes 
ambiguous. For example, if an X error occurs on both 



qubits one and two, then the syndrome result will be 1 1 ) . 
This will cause us to mis-correct by applying an X gate 
to qubit 3. Therefore, two errors will induce a logical bit 
flip and causes the code to fail, as expected. 

To be absolutely clear on how QEC acts to restore the 
system and protect against errors. Let us now consider 
a different and more physically realistic error mapping. 
We will assume that the errors acting on the qubits are 
coherent rotations of the form U = exjp(iea x ) on each 
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FIG. 3 Circuit required to encode and correct for a single X-error. We assume that after encoding a single bit-flip occurs on 
one of the three qubits (or no error occurs). Two initialized ancilla are then coupled to the data block which only checks the 
parity between qubits. These ancilla are then measured, with the measurement result indicating where (or if) an error has 
occurred, without directly measuring any of the data qubits. Using this syndrome information, the error can be corrected with 
a classically controlled X gate. 



qubit, with e <C 1. We choose coherent rotations so that 
we can remain in the state vector representation. This 
is not a necessary requirement , however more general in- 
coherent mappings would require us to move to density 
matrices. 

We assume that each qubit experiences the same error, 
hence the error operator acting on the state is, 

\il>) E = E\1>) L , 

E = U m = (cos(e)Z + i sin(e)cr x )® 3 

= colli + cx{o x <jioi + Gia x ai + <Ti<ti<t x ) 
+ Ci{o x a x ai + <7ia x a x + a x aicr x ) 
+ c 3 a x a x cr x . 

(21) 



where, 



c = cos 3 (e), 
c\ = icos 2 (e) sin(e) 
C2 = — cos(e) sin 2 (e) 
C3 = — zsin 3 (e). 



(22) 



Now let's examine the transformation that occurs when 
we run the error correction circuit in Fig. [3] which we 
denote via the unitary transformation, UqeCi over both 
the data and ancilla qubits, 



UqecE \i)) L |00) = c |V) £ |00) + c^ajo-j \^) L |11) + aja x ai \ijj) L |10) + *i<tio- x \ip) L |01>) 

+ c 2 {o- x o- x o- I \ip) L |01) + cri<J x a x \ip) L 1 1 1> + o x oio x \ip) L 1 10)) + c 3 <r x a x a x \ip) L |00) 



(23) 



Once again, the ancilla block is measured and the appro- 
priate correction operator is applied, yielding the results 



(up to renormalization) , 



Ancilla Measurement: 
Ancilla Measurement: 
Ancilla Measurement: 
Ancilla Measurement: 



1 00) , Collapsed State (with correction) 

1 01) , Collapsed State (with correction) 

1 10) , Collapsed State (with correction) 

1 11) , Collapsed State (with correction) 

I 



Co 



c 3 a x <r x a x \ip) L 
C20- X a x a x \ip) L 
C20- X a x a x \tp) L 
c%a x a x a x \tp) L 



(24) 



In each case, after correction (based on the syndrome result), we are left with approximately the same state. A 
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superposition of a "clean state" with the logically flipped 
state, o x a x a x \ip). The important thing to notice is the 
amplitudes related to the terms in the superposition. If 
we consider the unitary U acting on a single, uncncoded 
qubit, the rotation takes, 

U\i/>) = cos(e) | ip) + i sin(eK | j>) , (25) 

Consequently, the fidelity of the single qubit state is, 



F, 



unenco 



dcd = |(^|t/|^)| 2 = cos 2 



l-e 2 



(26) 



In contrast, the fidelity of the encoded qubit state after 
a cycle of error correction is, 



cos 6 (e) 



|c | 2 + |c 3 | 2 cos 6 (e)+sin 6 (e) (27) 



l-e 6 



with probability 1 — 3e 2 + 0(e 4 ) and 



ci 



error detected 



|ci| 2 +|c 2 | 2 



cos 4 (e) sin 2 (e) (28) 
cos 4 (e) sin 2 (e) + sin 4 (e) cos 2 (e) 
«l-e 2 . 

with probability 3e 2 +(9(e 4 ). This is the crux of how QEC 
suppresses errors at the logical level. During a round of 
error correction, if no error is detected (which if the er- 
ror rate is small, occurs with high probability) , the error 
on the resulting state is suppressed from 0(e 2 ) to 0(e 6 ), 
while if a single error is detected, the fidelity of the re- 
sulting state remains the same. 

This is expected, as the 3-qubit code is a single error 
correcting code. If one error has already been corrected 
then the failure rate of the logical system is conditional 
on experiencing one further error (which will be propor- 
tional to e 2 ). As e < 1 the majority of correction cycles 
will detect no error and the fidelity of the resulting en- 
coded state is higher than when unencoded. Note, that 
as e 2 — * 1/3 the benefit of the code disappears as every 
correction cycle detects an error and the resulting fidelity 
is no better than an unencoded qubit 

It should be stressed that no error correction 
scheme will, in general, restore a corrupted state 
to a perfectly clean code-state. The resulting state 
will contain a superposition of a clean state and cor- 
rupted states, the point is that the fidelity of the cor- 
rupted states, at the logical level, is greater than the cor- 
responding fidelity for unencoded qubits. Consequently 
the probability of measuring the correct result at the end 
of a specific algorithm increases when the system is en- 
coded. 

This example shows the basic principles of error cor- 
rection. As mentioned earlier, the 3-qubit code does not 
represent a full quantum code and the error model that 
we considered neglected imperfect gates and the possibil- 
ity of errors occurring during state preparation and/or 



correction. In the coming sections we will briefly take a 
look at several full quantum codes, both used for quan- 
tum memory and computation and we will introduce the 
concept of full QEC using stabilizer codes. This will then 
lead to a description of full fault-tolerant quantum error 
correction. 



V. THE NINE QUBIT CODE: THE FIRST FULL 
QUANTUM CODE 

The nine qubit error correcting code was first devel- 
oped by Shor (|Sho95|l in 1995 and is based largely on the 
3-qubit repetition code. The Shor code is a degenerate 
single error correcting code able to correct a logical qubit 
from one discrete bit flip, one discrete phase flip or one of 
each on any of the nine physical qubits and is therefore 
sufficient to correct for any continuous linear combina- 
tion of errors on a single qubit. The two basis states for 
the code are, 



1 
1 



\0) L = -^=(|000) + |111))(|000) + |111»(|000> + |111)) 



\1) L = ^=(1000) - |111))(|000> - |111»(|000) - |111)) 



and the circuit to perform the encoding is shown in Fig. [4] 
Correction for X errors, for each block of three qubits 
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FIG. 4 Circuit required to encode a single qubit with Shor's 
nine qubit code. 



encoded to (|000) ± |lll))/\/2 is identical to the three 
qubit code shown earlier. By performing the correction 
circuit shown in Fig.[3]for each block of three qubits, sin- 
gle a x = X errors can be detected and corrected. Phase 
errors (a z = Z) are corrected by examining the sign dif- 
ferences between the three blocks. The circuit shown 
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FIG. 5 Circuit required to perform phase correction for the 9-qubit code. 



in Fig. [5] achieves this. The first set of six CNOT gates 
compares the sign of blocks one and two of the code state 
and the second set of CNOT gates compares the sign for 
blocks two and three. Note that a phase flip on any one 
qubit in a block of three has the same effect, this is why 
the 9-qubit code is referred to as a degenerate code. In 
other error correcting codes, such as the 5- or 7-qubit 
codes (Stc96a; LMPZ96), there is a one-to-one mapping 
between correctable errors and unique states, in degener- 
ate codes such as this, the mapping is not unique. Hence 
provided we know in which block the error occurs it does 
not matter which qubit we apply the correction operator 
to. 

As the 9-qubit code can correct for single X errors in 
any one block of three and a single phase error on any 
of the nine qubits, this code is a full quantum error cor- 



recting code (we will detail in section IX why phase and 



bit correction is sufficient for the correction of arbitrary 
qubit errors) . Even if a bit and phase error occurs on the 
same qubit, the X correction circuit will detect and cor- 
rect for bit flips while the Z correction circuit will detect 
and correct for phase flips. As mentioned, the X error 
correction does have the ability to correct for up to three 
individual bit flips (provided each bit flip occurs in a dif- 
ferent block of three). However, in general the 9-qubit 
code is only a single error correcting code as it cannot 
handle multiple errors if they occur in certain locations. 

The 9-qubit code is in fact a member of a broader class 
of error correcting codes known as Bacon-Shor or subsys- 
tem codes ()Bac06[) . Subsystem codes have the property 



that certain subgroups of error operators do not corrupt 
the logical space. This can be seen by considering phase 
errors that occur in pairs for any block of three. For ex- 
ample, a phase flip on qubits one, two, four and five will 
leave both logical states unchanged. Subsystem codes are 
very nice codes from an architectural point of view. Error 
correction circuits and gates are generally simpler than 
for non-subsystem codes, allowing for circuit structures 
more amenable to the physical restrictions of a computer 
architecture (IAC07II . Additionally as subsystem codes 
that can correct for a larger number of errors have a simi- 
lar structure, we are able to perform dynamical switching 
between codes, in a fault-tolerant manner, which allows 
us to adapt the error protection in the computer to be 
changed depending on the noise present at a physical 
level (SEDH07). We will return and revisit subsystem 
codes later in section IXIV.AI 



VI. QUANTUM ERROR DETECTION 

So far we have focused on the ability to not only de- 
tect errors, but also to correct them. Another approach is 
to not enforce the correction requirement. Post-selected 
quantum computation, developed by Knill (Kni05) 
demonstrated that large scale quantum computing could 
be achieved with much higher noise rates when error de- 
tection is employed instead of more costly correction pro- 
tocols. The basic idea in post-selected schemes is to en- 
code the computer with error detecting circuits and if er- 
rors are detected, the relevant subroutine of the quantum 
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algorithm is reset and run again, instead of performing 
active correction. One of the downside to these types 
of schemes is that although they lead to large tolerable 
error rates, the resource requirements are unrealistically 
high. 

The simplest error detecting circuit is the 4-qubit 
code (1GBP97P . This encodes two logical qubits into four 
physical qubits with the ability to detect a single error 
on either of the two logical qubits. The four basis states 
for the code are, 



|00) = 


-^=(10000) - 


r 11111)), 


|01) = 


^(liioo)- 


-10011)), 


|10) = 


^(|ioio>- 


-|0101)), 


111) = 


^(lono)- 


- |1001)). 



(30) 



Fig. [6] illustrates the error detection circuit that can be 
utilized to detect a single bit and/or phase flip on one of 
these encoded qubits. If a single bit and/or phase flip oc- 
curs on one of the four qubits then the ancilla qubits will 
be measured in the |1) state. For example, let us consider 
the cases when a single bit flip occurs on one of each of 
the four qubits. The state of the system, just prior to the 
measurement of the ancilla is in table. In] Regardless of 



Error Location 


Final State, |data) | ancilla) 


No Error 


WJOO) 


Qubit 1 


Xx |V) £ 


10) 


Qubit 2 




10) 


Qubit 3 




10) 


Qubit 4 




10) 



TABLE II Qubit and ancilla state, just prior to measurement 
for the 4-qubit error detection code when a single bit-flip has 
occurred on at most one of the four qubits. 



the location of the bit flip, the ancilla system is measured 
in the state |10). Similarly if one considers a single phase 
error on any of the four qubits the ancilla measurement 
will return |01). In both cases no information is obtained 
regarding where the error has occurred, hence it is not 
possible to correct the state. Instead the subroutine can 
be reset and re-run. 



VII. STABILIZER FORMALISM 

So far we have presented error correcting codes from 
the perspective of their state representations and their 
preparation and correction circuits. This is a rather in- 
efficient method for describing the codes as the state 
representations and circuits clearly differ from code to 
code. The majority of error correcting codes that are 
used within the literature are members of a class known 



as stabilizer codes. Stabilizer codes are very useful to 
work with. The general formalism applies broadly and 
there exists general rules to construct preparation cir- 
cuits, correction circuits and fault-tolerant logical gate 
operations once the stabilizer structure of the code is 
specified. 

The stabilizer formalism which was first introduced by 
Daniel Gottesman (|Got97jl uses essentially the Heisen- 
berg representation for quantum mechanics which de- 
scribes quantum states in terms of operators rather that 
the basis states themselves. An arbitrary state \tp) is de- 
fined to be stabilized by some operator, K, if it is a +1 
eigenstate of K, i.e. 



K\r/>) = \il>). 

For example, the single qubit state 
the operator K = er z , i.e. 



(31) 

is stabilized by 
(32) 



Defining multi-qubit states with respect to this formalism 
relics on the group structure of multi-qubit operators. 

Within the group of all possible, single qubit operators, 
there exists a subgoup, denoted the Pauli group, V, which 
contains the following elements, 

V = {±a I ,±ia I ,±a x ,±ia x ,±cr y ,±ia y ,±a z ,±ia z }. 

(33) 

It is easy to check that these matrices form a group un- 
der multiplication through the commutation and anti- 
commutation rules for the Pauli set, {<Ji\ — {a x ,a y ,a z }, 



[<Ti,(Tj] = 2ieijk(Tk, {cj, (Tj} = 2di 



(34) 



where, 



and 



+lfor fc)e {(1,2,3), (2, 3,1), (3, 1,2)} 
-lfor fc)e {(1,3, 2), (3, 2,1), (2, 1,3)} 
for i = j, j = k, or k = i 

(35) 



I 1 for i = j 
° 11 ~) for i ^ j. 



(36) 



The Pauli group extends over N-qubits by simply taking 
the N fold tensor product of V, i.e. 

V N = V® N 

= {±cr/, ±i(jj, ±a x , ±ia x , ±a y , ±i<7 y ,±a z , ±ia z }® N . 

(37) 

An iV-qubit stabilizer state, \i/j) n is then defined via an 
A^-element Abelian subgroup, Q, of the A^-qubit Pauli 
group, in which \ip) N is a +1 eigenstate of each element, 



{Gi\Gi M = |V) , [Gi, Gj] = V c V N . 



(38) 
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FIG. 6 Circuit required detect errors in the 4-qubit error detection code. If both ancilla measurements return |0), then the 
code state is error free. If either measurement returns |1), an error has occurred. Unlike the 9-qubit code, the detection of an 
error does not give sufficient information to correct the state. 



Given this definition, the state \ip) N can be equivalcntly 
defined either through the state vector representation or 
by specifying the stabilizer set, Q. 

Many extremely useful multi-qubit states are stabi- 
lizer states, including two-qubit Bell states, Greenberger- 
Horne-Zeil inger (GHZ) s tates (|GHZ89I IGHSZ90|) . Clus- 
ter states (IBROli IRB01|) and codeword states for QEC. 
As an example, consider a three qubit GHZ state, defined 
as, 



IGHZ) 



|000) 



1111) 



V2 



(39) 



This state can be expressed via any three linearly in- 
dependent elements of the |GHZ) 3 stabilizer group for 
example, 



G\ = <t x ® <r x ® a x = XXX, 

G 2 = CT z (g)<T z (g)(J I = ZZI, 

G 3 = <7j ® <j z ® a z = IZZ. 



(40) 



where the right-hand side of each equation is the short- 
hand representation of stabilizers. Note that these three 
operators form an Abelian group [Eq. 38 as, 



W = G t G 3 \i>) - 

= m - m 



GjGi |V) 
= 0, V 



(41) 



Similarly, the four orthogonal Bell states, 



±\_ |00)±|11) 



I* > = 



±x _ |01)±|10) 



V2 
V2 



(42) 



are stabilized by the operators, G\ = (—l) a XX, and 
G 2 = (-l) b ZZ. Where [a, b] € {0,1} and each of the 
four Bell states correspond to the four unique pairs, 
{$+,*+,*-,*-} = {[0,0], [0,1], [1,0], [1,1]}. 
VIII. QEC WITH STABILIZER CODES 

The use of the stabilizer formalism to describe quan- 
tum error correction codes is extremely useful since it 
allows for easy synthesis of correction circuits and also 
clearly shows how logical operations can be performed 
directly on encoded data. As an introduction we will fo- 
cus on arguably the most well known quantum code, the 
7-qubit Steane code, first proposed in 1996 ()Ste96ajl . 

The 7-qubit code represents a full quantum code that 
encodes seven physical qubits into one logical qubit, with 
the ability to correct for a single X and/or Z error. The 
|0) L and |1) L basis states are denned as, 



|0) L = ^(|0000000) + |1010101) + 10110011) + 11100110) + 10001111) + 11011010) + 10111100) + 11101001)), 

\1) L = -pfllllllll) + |0101010) + 11001100) + 10011001) + 11110000) + |0100101) + |1000011) + 10010110)). 
V8 



(43) 
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The stabilizer set for the 7-qubit code is fully specified 
by the six operators, 



K 1 = IIIXXXX, 
K 3 = IXXIIXX, 



K z = XIXIXIX, 

K A = IIIZZZZ (44) 



K 5 = ZIZIZIZ, K b = IZZIIZZ. 

As the 7-qubit codeword states are specified by only 
six stabilizers, the code contains two basis states, which 
are the logical states. With a final operator, K 7 — 
ZZZZZZZ = Z® 7 fixing the state to one of the code- 
words, K 7 \Q) L = |0) L and K 7 \l) L = -\1) L . The 7- 
qubit code is defined as a [[n, k,d]] = [[7,1,3]] quantum 
code, where n = 7 physical qubits encode k = 1 logical 
qubit with a distance between basis states d = 3, correct- 
ing t = (d — l)/2 = 1 error. Notice that the stabilizer set 
separates into X and Z sectors which defines the code 
as a Calderbank-Shor-Steane (CSS) code. CSS codes are 
extreamly useful since they allow for straightforward log- 
ical gate operations to be applied directly to the encoded 
data [Section XI and are reasonably easy to derive from 
classical codes. 

Although the 7-qubit code is the most well known 
Stabilizer code, there are other stabilizer codes which 
encode multiple logical qubits and correct for more er- 
rors (|Got97p . The downside to these lager codes is that 
they require more physical qubits and more complicated 
error correction circuits. Tables IIIII and IIVI shows the 
stabilizer structure of two other well known codes, the 
9-qubit code (lSho95i) which we have examined and the 
5-qubit code (LMPZ96) which represents the smallest 
possible quantum code that corrects for a single error. 
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X 
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TABLE III The eight Stabilizers for the 9-qubit Shor code, 
encoding nine physical qubits into one logical qubit to correct 
for a single X and/or Z error. 



A. State Preparation 

Using the stabilizer structure for QEC codes, the log- 
ical state preparation and error correcting procedure is 
straightforward. Recall that the codeword states are de- 
fined as +1 eigenstates of the stabilizer set. In order 
to prepare a logical state from some arbitrary input, we 
need to forcibly project qubits into eigenstates of these 
operators. 
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TABLE IV The Four Stabilizers for the [[5,1,3]] quantum 
code, encoding five physical qubits into one logical qubit to 
correct for a single X and/or Z error. Unlike the 7- and 9- 
qubit codes, the [[5,1,3]] code is a non-CSS code, since the 
stabilizer set does not separate into X and Z sectors. 



Consider the circuit shown in Fig. [7] For some arbi- 
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FIG. 7 Quantum Circuit required to project an arbitrary 
state, \ip) j into a ±1 eigenstate of the Hermitian operator, 
U = . The measurement result of the ancilla determines 
which eigenstate is projected to. 



trary input state, \ip)j, an ancilla which is initialized in 
the |0) state is used as a control qubit for a Hermitian 
operation (W = U) on \ip)r. After the second Hadamard 
gate is performed, the state of the system is, 



1 



(|V>j + 0»j) 



j(|V>i-l/W>>/)|l>. (45) 



The ancilla qubit is then measured in the computational 
basis. If the result is |0), the input state is projected to 
(neglecting normalization), 



\^) F = \i>) I + u\ 



1 1 ■ 



(46) 
is a +1 



Since U is Hermitian, U F = \ip) F , hence \ip) F 
eigenstate of U. If the ancilla is measured to be |1), then 
the input is projected to the state, 



\i>) F = \i>)j-u\ 



1 1 ' 



(47) 



which is the —1 eigenstate of U. Therefore, provided U 
is Hermitian, the general circuit of Fig. [7] will project an 
arbitrary input state to a ±1 eigenstate of U. This pro- 
cedure is well known and is refered to as either a "parity" 
or "operator" measurement (jNCOOp . 

From this construction it should be clear how QEC 
state preparation proceeds. Taking the [[7, 1, 3]] code 
as an example, 7-qubits are first initialized in the state 
|0)® 7 , after which the circuit shown in Fig. [t] is applied 
three times with U = (K 1 , K 2 ,K 3 ), projecting the input 
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state into a simultaneous ±1 eigenstate of each X stabi- 
lizer describing the [[7,1,3]] code. The result of each 
operator measurement is then used to classically con- 
trol a single qubit Z gate which is applied to one of the 
seven qubits at the end of the preparation. This single Z 
gate converts any —1 projected eigenstates into +1 eigen- 
states. Notice that the final three stabilizers do not need 
to be measured due to the input state, |0)® 7 , already be- 
ing a +1 eigenstate of (if 4 , K 5 , if 6 ). Fig. 8] illustrates 
the final circuit, where instead of one ancilla, three are 
utilized to speed up the state preparation by performing 
each operator measurement in parallel. 
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FIG. 8 Quantum circuit to prepare the [[7,1,3]] logical |0) 
state. The input state |0}® 7 is projected into an eigenstate of 
each of the X stabilizers shown in Eq. |44| After each ancilla 
measurement the classical results are used to apply a single 
qubit Z gate to qubit i = l Mz +2 Ms +4 Ml which converts the 
state from a —1 eigenstates of (K , K 2 , if 3 ) to +1 eigenstates. 



As a quick aside, let us detail exactly how the rela- 
vant logical basis states can be derived from the stabilizer 
structure of the code by utilizing the preparation circuit 
illustrated above. Instead of the 7-qubit code, we will 
use the stabilizer set shown in Table IIVI to calculate the 
|0)r state for the 5-qubit code. The four code stabilizers 
are given by, 



if 1 = XZZXI, 
if 3 = XIX ZZ, 



if 2 = IXZZX, 
4 ( 48 ) 



K 



ZXIXZ. 



As with the 7-qubit code, projecting an arbitrary state 
into a 4-1 eigenstate of these operators define the two, 
logical basis states |0) L and |1) L , with the operator 
Z = ZZZZZ, fixing the state to either |0) L or \l) L . 
Therefore, calculating |0) L from some initial un-encoded 
state requires us to project the initial state into a 4-1 
eigenstate of these operators. If we take the initial, un- 
encoded state as |00000), then it already is a 4-1 eigen- 
state of Z. Therefore, to find |0) L we simply calculate, 



|0) L = ]J(/® 5 + ir) 100000) 



(49) 



up to normalization, 
find, 

1 



Expanding out this product, we 



( 100000) 4- 101010) 4- |10100) - |11110) 4- 

|01001) - |00011) - |11101) - |10111) + 
|10010) - |11000) - |00110) - |01100) - 
|11011) - |10001) - |01111) 4- |00101)). 

(50) 

Note, that the above state vector does not match up with 
those given in (L MPZ96 ,) . However, these vectors are 
equivalent up to local rotations on each qubit. Therefore, 
matching up the original state requires locally perturbing 
the stabilizer set to reflect these rotations. 



B. Error Correction 

Error correction using stabilizer codes is a straightfor- 
ward extension of state preparation. Consider an arbi- 
trary single qubit state that has been encoded, 



a|0) + /3|l) 



(51) 



Now assume that an error occurs on one (or multiple) 
qubits which is described via the operator E, where E is 
a combination of X and/or Z errors over the N physi- 
cal qubits of the logical state. By definition of stabilizer 
codes, K % \i>) L = \ip) L , i € [1, ..,N — k], for a code en- 
coding k logical qubits. Hence the erred state, E 
satisfies, 



L ' 



K l E\ 



{-vrEK i \^) L = {-vrE\^) 1 



(52) 



i=l 



where m is defined as, m — 0, if [E, if 1 ] = and m = 1, 
if {E, K 1 } — 0. Therefore, if the error operator com- 
mutes with the stabilizer, the state remains a +1 eigen- 
state of K l , if the error operator anti-commutes with the 
stabilizer then the logical state is flips to now be a — 1 
eigenstate of K % . 

Hence the general procedure for error correction is 
identical to state preparation. Each of the code stabi- 
lizers are sequentially measured. Since a error free state 
is already a 4-1 eigenstate of all the stabilizers, any er- 
ror which anti-commutes with a stabilizer will flip the 
eigenstate and consequently the parity measurement will 
return a result of |1). 

Taking the [[7,1,3]] code as an example, you can see 
that if the error operator is E = Xi, where i = (1,...,7), 
representing a bit-flip on any one of the 7 physical qubits, 
then regardless of the location, E will anti-commute with 
a unique combination of (K 4 , K 5 , K 6 ). Hence the classi- 
cal results of measuring these three operators will indi- 
cate if and where a single X error has occurred. Similarly, 
if E = Zi, then the error operator will anti-commute with 
a unique combination of, (iT 1 , if 2 , if 3 ). Consequently, 
the first three stabilizers for the [[7, 1, 3]] code correspond 
to Z sector correction while the second three stabilizers 
correspond to X sector correction. Note, that correction 
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for Pauli Y errors are also taken care of by correcting in 
the X and Z sector since a Y error on a single qubit is 
equivalent to both an X and Z error on the same qubit, 
i.e. Y = iXZ. 

Fig. [9] illustrates the circuit for full error correction 
with the [[7,1,3]] code. As you can see it is simply an 
extension of the preparation circuit [Fig. [8] where all 
six stabilizers are measured across the data block. Even 
though we have specifically used the [[7,1,3]] code as 
an example, the procedure for error correction and state 
preparation is identical for all stabilizer codes allowing 
for full correction for both bit and phase errors without 
obtaining any information regarding the state of the log- 
ical qubit. 



IX. DIGITIZATION OF QUANTUM ERRORS 

Up until now we have remained fairly abstract regard- 
ing the analysis of quantum errors. Specifically, we have 
examined QEC from the standpoint of a discrete set of 
Pauli errors occurring at certain locations within a larger 
quantum circuit. In this section we examine how this 
analysis of errors relates to more realistic processes such 
as environmental decoherence and systematic gate errors. 

Digitization of quantum noise is often assumed when 
people examine the stability of quantum circuit design 
or attempt to calculate thresholds for concatenated error 
correction. However, the equivalence of discrete Pauli er- 
rors to more general, continuous, noise only makes sense 
when we consider the stabilizer nature of the correction 
procedure. Recall from section|VTT]that correction is per- 
formed by re-projecting a potentially corrupt data block 
into +1 eigenstates of the stabilizer set. It is this pro- 
cess that acts to digitize quantum noise, since a general 
continuous mapping from a "clean" codeword state to 
a corrupt one will not satisfy the stabilizer conditions, 
we will first introduce how a coherent systematic error, 
caused by imperfect implementation of quantum gates, 
are digitized during correction, after which we will briefly 
discuss environmental decoherence from the standpoint 
of the Markovian decoherence model. 



where \ip') L = Um |^)r is the perfectly applied N qubit 
gate, (i.e. the stabilizer set for T remains invariant 
under the operation Um [see Sec. XI ). and Ue is a co- 



herent error operator which is expanded in terms of the 
N qubit Pauli Group, Ej <G Pm- Now append two ancilla 

blocks, |^4o) X and |^4o) Z j which are all initialized and are 
used for X and Z sector correction, then run a full error 
correction cycle, which we represent by the unitary oper- 
ator, Uqec- It will be assumed that \i/>) L is encoded with 
a hypothetical QEC code which can correct for N errors 
(both X and/or Z), hence there is a one-to-one mapping 
between the error operators, Ej, and the orthogonal basis 
states of the ancilla blocks, 



\A ) X \A ) Z 



UqecU'm I V/ jj 1 

= U QECz 2 a J E iW)L\Ao) X \A Q ) 

3 

= Y,«i E iW)L\Ai) x \Ai) z - 



(54) 



The ancilla blocks are then measured, projecting the data 
blocks into the state Ej \ip') L with probability \ctj\ 2 , after 
which the correction Ej is applied based on the syndrome 
result. As the error operation Ej is simply an element of 
Vn, correcting for X and Z independently is sufficient to 
correct for all error operators (as Y errors are corrected 
when a bit and phase error is detected and corrected on 
the same qubit). 

For very small systematic inaccuracies, the expansion 
co-efficient, ao, which corresponds to E Q — I® N will be 
very close to 1, with all other co-efficients small. Hence 
during correction there will be a very high probability 
that no error is detected. This is the digitization effect of 
quantum error correction. Since codeword states are spe- 
cific eigenstates of the stabilizers, then the re-projection 
of the state when each stabilizer is measured forces any 
continuous noise operator to collapse to the discrete Pauli 
set, with the magnitude of the error dictating the prob- 
ability that the data block is projected into discrete per- 
turbation of a "clean" state. 



A. Systematic gate errors 

We have already shown an example of how systematic 
gate errors are digitized into a discrete set of Pauli opera- 
tors in Sec. |III| However, in that case we only considered 
a very restrictive type of error, namely the coherent oper- 
ator U = cxp(ieX). We can easily extend this analysis to 
cover all forms of systematic gate errors. Consider an N 
qubit unitary operation, Um, which is valid on encoded 
data. Assume that Um is applied inaccurately such that 
the resultant operation is actually U' N . Given a general 
encoded state \i/>) N , the final state can be expressed as, 

U' N W) l = UeUm \4>) l = ]T ocjEj \il/) L , (53) 

3 



B. Environmental decoherence 

A complete analysis of environmental decoherence in 
relation to quantum information is a lengthy topic. In- 
stead of a detailed review, we will instead simply present 
a specific example to highlight how QEC relates to envi- 
ronmental effects. 

The Lindblad formalism (jGar911 INC00I IDWM03[) pro- 
vides an elegant method for analyzing the effect of deco- 
herence on open quantum systems. This model does have 
several assumptions, most notably that the environmen- 
tal bath couples weakly to the system (Born approxima- 
tion) and that each qubit experiences un-correlated noise 
(Markovian approximation). While these assumptions 
are utilized for a variety of systems (|BHPC03I BM03 
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FIG. 9 Quantum circuit to to correct for a single X and/or Z error using the [[7, 1,3]] code. Each of the 6 stabilizers are 
measured, with the first three detecting and correcting for Z errors, while the last three detect and correct for X errors. 



BKDQ4]), it is k nown that t hey may not hold in some 
cases (|HMCS00l IMCM+05I IAPN+05I IALKH02p . Par- 
ticularly in superconducting systems where decoherence 
can be caused by small numbers of fluctuating charges. 
In this case more specific decoherence models need to be 
considered. 

Using this formalism, the evolution of the density ma- 
trix can be written as, 

d tP ^- l -[H,p} + J2^kC[p}. (55) 

k 

Where H is the Hamiltonian, representing coherent, dy- 
namical evolution of the system and Ck[p] = {[Lk, pLU + 
[LkP, l\.])/2 represents the incoherent evolution. The op- 
erators Lfc are known as the Lindblad quantum jump 
operators and are used to model specific decoherence 
channels, with each operator parameterized by some rate 
Tfc > 0. This differential equation is known as the quan- 
tum louiville equation or more generally, the density ma- 
trix master equation. 

To link Markovian decoherence to QEC, consider a 
special set of decoherence channels that help to sim- 
plify the calculation, representing a single qubit un- 
dergoing dephasing, spontaneous emission and sponta- 
neous absorption. Dephasing of a single qubit is mod- 
elled by the Lindblad operator L\ = Z while sponta- 
neous emission/ absorption are modelled by the operators 
L 2 = |0) (1| and L 3 — |1) (0| respectively. For the sake 
of simplicity we assume that absorption/emission occur 
at the same rate, T. Consequently, the density matrix 
evolution is given by, 



d t p=~-[H,p]+T z (ZpZ-p) 



; (XpX + YpY-2p). 

(56) 

If it is assumed that the qubit is not undergoing any 
coherent evolution (H = 0), i.e. a memory stage within 



a quantum algorithm, then Eq. [56] can be solved by re- 
expressing the density matrix in the Bloch formalism. 
Set p(t) = 1/2 + x(t)X + y(t)Y + z(t)Z, then Eq. [56] 
with H = 0, reduces to, d t S(t) = AS(t) with S(t) = 
(x(t),y(t),z(t)) T and 



.4 




(57) 



This differential equation is easy to solve, leading to, 



p(t) = [l-p(t)}p(0) +Px (t)Xp(0)X 
+ p y (t)Yp(0)Y + Pz (t)Zp(0)Z, 



(58) 



where, 



1 



Px(t) =Py(t) = j(l 
1 



-2Tt 



-in _ 2 e -( r + 2r *) t ) j 

P(t)=p x (t)+Py{t)+ Pz {t). 



p z (t)=-(l + e- 



(59) 



If this single qubit is part of a QEC encoded data block, 
then each term represents a single error on the qubit ex- 
periencing decoherence. Two blocks of initialized ancilla 
qubits are added to the system and the error correction 
protocol run. Once the ancilla qubits are measured, the 
state will collapse to no error, with probability 1 — p(t), 
or a single X,Y or Z error, with probabilities p x (t),p y (t) 
and p z (t). 

We can also see how temporal effects are incorporated 
into the error correction model. The temporal integra- 
tion window t of the master equation will influence how 
probable an error is detected and corrected for a fixed 
rate T. The longer between correction cycles, the more 
probable the qubit experiences an error. 



17 



C. More General mappings 

Both the systematic gate errors and the errors induced 
by environmental decoherence illustrate the digitization 
effect of quantum error correction. However, we can quite 
easily generalize digitization to arbitrary mappings of the 
density matrix. In this case consider a more general 
Krauss map on a multi-qubit density matrix, 

p^J2 A kP A * (60) 

k 

where ^A^Ak = !■ For the sake of simplicity let us 
choose a simple mapping where A-y = (Zi + iZ 2 ) / 'y/2 and 
Ak = for k ^ 1. This mapping essentially represents 
dephasing on two qubits. However, this type of map- 
ping (when considered in the context of error correction) 
represents independent Z errors on cither qubit one or 



two. 

To illustrate, first expand out the density matrix (ne- 
glecting normalization) , 

p — * A\pAi = Z\pZ\ + Z 2 pZ 2 - iZipZ 2 + iZ 2 pZx (61) 

Note that only the first two terms in this expansion, on 
their own, represent physical mixtures, the last two off- 
diagonal terms are actually irrelevant in the context of 
QEC and are removed during correction. To illustrate we 
again assume that p represents a protected qubit, where 
Z\ and Z 2 are physical errors on qubits comprising the 
codeblock. As we are only considering phase errors in this 
example, we will ignore X correction (but the analysis 
automatically generalizes if the error mapping contains X 
terms). A fresh ancilla block, represented by the density 
matrix /?§ is coupled to the system and the unitary Uqec 
is run, 



u qecP' ® PoUqec =U i QEC Z 1 pZ 1 <g) pIUqec + U^ QEC Z 2 pZ 2 ® p z U QE c 

-iU f QEC Z lP Z 2 <8> p z U QEC + iUl }EC Z 2 pZ 1 ® plU QEC (62) 
=Z lP Z x <gs \Z X ) (Zi\ + Z 2 pZ 2 <g> \Z 2 ) (Z 2 \ - iZ lP Z 2 ® \Zi) (Z 2 \ + tZ 2P Z ± \Z 2 ) (Zi\ 



where \Z\) and \Z 2 ) represent the two orthogonal syn- 
drome states of the ancilla that are used to detect phase 
errors on qubits one and two respectively. The important 
part of the above expression is that when the syndrome 
qubits arc measured we are calculating Tr(p \Zi) {Zi\) or 
Tr(p \Z 2 ) (Z 2 \), therefore the two cross terms in the above 
expression are never observed. In this mapping the only 
two possible states that exist after the measurement of 
the ancilla system are, 

Z x pZ\ ® \Z\) {Zi\ with Probability = - 

\ (63) 

Z 2 pZ 2 <g) \ Z 2 ) (Z 2 \ with Probability = - 



X. FAULT-TOLERANT QUANTUM ERROR 
CORRECTION AND THE THRESHOLD THEOREM. 



Therefore, not only are the cross terms eliminated via er- 
ror correction but the final density matrix again collapses 
to a single error perturbation of "clean" codeword states 
with no correlated errors. 

Consequently, in standard QEC analysis it is assumed 
that after each elementary gate operation, measurement, 
initialization and memory step, a hypothetical error cor- 
rection cycle is run. This cycle digitizes all continuous 
errors (either systematic or environmental) into either 
an X and/or Z error on each qubit. This cycle is as- 
sumed to be error free and take zero time. In this way 
error correction can be analyzed by assuming perfect gate 
operations and discrete, probabilistic errors. The prob- 
ability of each error occuring can then be independently 
calculated via a systematic gate analysis or through the 
evolution of the master equation. 



Section VII detailed the protocols required to correct 
for quantum errors, however this implementation of QEC 
assumed the following, 

1. Errors only occur during "memory" regions, i.e. 
when quantum operations or error correction are 
not being performed and we assume errors do not 
occur on ancilla qubits. 

2. The quantum gates themselves do not induce any 
systematic errors within the logical data block. 

Clearly these are two very unrealistic assumptions and 
error correction procedures and logical gate operations 
need to be designed such that they can still correct for 
errors. 



A. Fault-tolerance 

The concept of Fault-tolerance in computation is not 
a new idea, it was firs t dev eloped in relation to classi- 
cal computing (|Neu55| IG83 Avi87). However, in recent 
years the precise manufacturing of digital circuitry has 
made large scale error correction and fault-tolerant cir- 
cuits largely unnecessary. 

The basic principle of Fault-tolerance is that the cir- 
cuits used for gate operations and error correction pro- 
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cedures should not cause errors to cascade. This can be 
seen clearly when we look at a simple CNOT operation 
between two qubits [Fig. 10 . In this circuit we are per- 



forming a sequence of three CNOT gates which act to 
take the state | 111) |000) -> 1 11 1) 1 11 1) . In Fig. Wok. we 



consider a single X error which occurs on the top most 
qubit prior to the first CNOT. This single error will cas- 
cade through each of the three gates such that the X error 
has now propagated to four qubits. Fig. [Top , shows a 
slightly modified design that implements the same oper- 
ation, but the single X error now only propagates to two 
of the six qubits. If we consider each block of three as a 
single logical qubit, then the staggered circuit will only 
induce a total of one error in each logical block, given a 
single X error occurred somewhere during the gate op- 
erations. This is the one of the standard definitions of 
Fault-tolerance. 

fault-tolerant circuit element: A single error will cause 
at most one error in the output for each logical qubit 
block. 

It should be stressed that the idea of Fault-tolerance 
is a discrete definition, either a certain quantum opera- 
tion is fault-tolerant or it is not. What is defined to be 
fault-tolerant can change depending on the error correc- 
tion code used. For example, for a single error correcting 
code, the above definition is the only one available (since 
any more than one error in a logical qubit will result in 
the error correction procedure failing). However, if the 
quantum code employed is able to correct multiple errors, 
then the definition of Fault-tolernace can be relaxed, i.e. 
if the code can correct three errors then circuits may 
be designed such that a single failure results in at most 
two errors in the output (which is then correctable). In 
general, for an code correcting t = [(d — 1)/2J errors, 
fault-tolerance requires that < t errors during an opera- 
tion does not result in > t errors in the output for each 
logical qubit. 



B. Threshold Theorem 

The threshold theorem is truly a remarkable result 
in quantum information and is a consequence of fault- 
tolerant circuit design and the ability to perform dy- 
namical error correction. Rather than present a detailed 
derivation of the theorem for a variety of noise models, 
we will instead take a very simple case where we utilize 
a quantum code that can only correct for a single er- 
ror, using a model that assumes uncorrelated, errors on 
individual qubits. For more rigorous derivations of the 
theorem see (IABQ97I l( lot 971 lAW]) . 

Consider a quantum computer where each physical 
qubit experiences either an X and/or Z error indepen- 
dently with probability p, per gate operation. Further- 
more, it is assumed that each logical gate operation and 
error correction circuit is designed according to the rules 
of Fault-tolerance and that a cycle of error correction is 
performed after each elementary logical gate operation. 



If an error occurs during a logical gate operation, then 
Fault-tolerance ensures this error will only propagate to 
at most one error in each block, after which a cycle of 
error correction will remove the error. Hence if the fail- 
ure probability of un-encoded qubits per time step is p, 
then a single level of error correction will ensure that the 
logical step fails only when two (or more) errors occur. 
Hence the failure rate of each logical operation, to lead- 
ing order, is now p\ = cp 2 , where p\ is the failure rate 
(per logical gate operation) of a 1st level logical qubit and 
c is the upper bound for the number of possible 2-error 
combinations which can occur at a physical level within 
the circuit consisting of the correction cycle + gate op- 
eration + correction cycle (Ali07). We now repeat the 
process, re-encoding the computer such that a level-2 
logical qubit is formed, using the same [[n, k, d]] quan- 
tum code, from n, level-1 encoded qubits. It is assumed 
that all error correcting procedures and gate operations 
at the 2nd level are self-similar to the level-1 operations 
(i.e. the circuit structures for the level-2 encoding are 
identical to the level- 1 encoding). Therefore, if the level- 
1 failure rate per logical time step is p\, then by the 
same argument, the failure rate of a 2-level operation is 
given by, p\ — c{p 1 L ) 2 — c 3 p 4 . This iterative procedure 
is then repeated (referred to as concatenation) up to the 
fcth level, such that the logical failure rate, per time step, 
of a fc-level encoded qubit is given by, 



(64) 



Eq. [64] implies that for a finite physical error rate, p, 
per qubit, per time step, the failure rate of the fcth-level 
encoded qubit can be made arbitrarily small by simply 
increasing k, dependent on cp < 1. This inequality de- 
fines the threshold. The physical error rate experienced 
by each qubit per time step must be pth < 1/c to ensure 
that multiple levels of error correction reduce the failure 
rate of logical components. 

Hence, provided sufficient resources are available, an 
arbitrarily large quantum circuit can be successfully im- 
plemented, to arbitrary accuracy, once the physical er- 
ror rate is below threshold. The calculation of thresh- 
olds is therefore an extremely important aspect to quan- 
tum architecture design. Initial estimates at the thresh- 
old, which gave p th ~ 10" 4 (!Kit97l IABQ97I IGot97l) did 
not sufficiently model physical syst ems in an accurate 
way. Rece nt results (|SFH08l ISDT07I ISBF+06I IMCT+04I 
BKSO05) have been estimated for more realistic quan- 
tum processor architectures, showing significant differ- 
ences in threshold when architectural considerations are 
taken into account. 



XI. FAULT-TOLERANT OPERATIONS ON ENCODED 
DATA 



Sections | VII| and |X] showed how fault-tolerant QEC al- 
lows for any quantum algorithm to be run to arbitrary 
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FIG. 10 Two circuits to implement the transformation | 111) |000) — ♦ | 111) | 111) . a) shows a version where a single X error can 
cascade into four errors while b) shows an equivalent circuit where the error only propagates to a second qubit. 



accuracy. However, the results of the threshold theorem 
assume that logical operations can be performed directly 
on the encoded data without the need for continual de- 
coding and re-encoding. Using stabilizer codes, a large 
class of operations can be performed on logical data in 
an inherently fault-tolerant way. 

If a given logical state, \ip) L , is stabilized by K, and 
the logical operation U is applied, the new state, U \i/)) L 
is stabilized by UKU\ i.e, 



UKU^U = UK \ip) L = U 



(65) 



In order for the codeword states to remain valid, the sta- 
bilizer set for the code, Q, must remain fixed through 
every operation. Hence for U to be a valid operation on 
the data, UQW = Q. 



A. Single Qubit Operations 

The logical X and Z operations on a single encoded 
qubit are the first examples of valid codeword operations. 
Taking the [[7,1,3]] code as an example, X and Z are 
given by, 

X = XXXXXXX = X® 7 , z = zzzzzzz = z® 7 . 

(66) 

Since the single qubit Pauli operators satisfy XZX — —Z 
and ZXZ = -X then, XK l X = K { and ZK l Z = K l for 
each of the [[7, 1, 3]] stabilizers given in Eq. 
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The fact 

that each stabilizer has a weight of four guarantees that 
UKU^ picks up an even number of —1 factors. Since the 
stabilizers remain fixed the operations are valid. How- 



ever, what transformations do Eq. 66 actually perform 
on encoded data? 

For a single qubit, a bit-flip operation X takes |0) <-> 
|1). Recall that for a single qubit Z |0) = |0) and 



Z\l) = — |1), hence for X to actually induce a logi- 



cal bit- flip it must take, |0) i 



1 



For the [[7,1,3]] 



code, the final operator which fixes the logical state is 



K 7 = Z® 7 , where K 7 \Q) 1 



\Q) L and ^ 7 |1) L 



1 



As XK 7 X = —K 7 , any state stabilized by K 7 becomes 
stabilized by — K 7 (and vice-versa) after the operation of 
X. Therefore, X represents a logical bit flip. The same 
argument can be used for Z by considering the stabilizer 
properties of the states |±) = (|0) ± |l))/v2- Hence, the 
logical bit- and phase-flip gates can be applied directly 
to logical data by simply using seven single qubit X or 
Z gates, [Fig. 
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Two other useful gates which can be applied in this 
manner is the Hadamard rotation and phase gate, 



H = 



1 

71 



P = 



1 

i 



(67) 



These gates are useful since when combined with the 
two-qubit CNOT gate, they can generate a subgroup 
of all multi-qubit gates known as the Clifford group 
(gates which map Pauli group operators back to the Pauli 
group). Again, using the stabilizers of the [[7, 1, 3]] code 
and the fact that for single qubits, 



HXH = Z, 
PXP^ = iXZ, 



HZH = X, 
PZP^ = Z, 



(68) 



a seven qubit bit-wise Hadamard gate will switch X 
with Z and therefore will simply flip {K 1 , K 2 , K 3 } with 
{K 4 , K 5 , K 6 }, and is a valid operation. The bit-wise 
application of the P gate will leave any Z stabilizer in- 
variant, but takes X — > iXZ. This is still valid since 
provided there are a multiple of four non-identity oper- 
ators for the stabilizer set, the factors of i will cancel. 
Hence seven bit- wise P gates is valid for the [[7, 1, 3]] 
code. 
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FIG. 11 Bit-wise application of single qubit gates in the [[7,1,3]] code. Logical X, Z H and P gates can trivially be applied 
by using seven single qubit gates, fault-tolerantly. Note that the application of seven P gates results in the logical P+ being 
applied and vice-versa. 



What does H and P do to the logical state? For a sin- 
gle qubit, the Hadamard gate flips any Z stabilized state 
to a X stabilized state, i.e |0, 1) <-» |+, — ). Looking at 
the transformation of K 7 , HK 7 H — X® 7 , the bit-wise 
Hadamard gate will invoke a logical Hadamard operation. 
The single qubit P gate leaves a Z stabilized state invari- 
ant, while an X eigenstate becomes stabilized by iXZ. 
Hence, Pt(X® 7 )P = -i(XZ)® 7 and the bit-wise gate, P, 
represent^ a logical P^ gate on the data block. Similarly, 

Each 



11 



bit-wise Pt gates enact a logical P gate [Fig. 
of these fault-tolerant operations on a logically encoded 
block are commonly referred to as transversal operations, 
as a logical operation is obtained by a set of individual 
operations acting transversally on the physical qubits. 



the two qubit stabilizer set, 

X <g> I -> X <g> X, 
I (g) Z -> Z <g) Z, 
Z®I->Z®I, 
I®X^I®X. 



(69) 



Where the first operator corresponds to the control qubit 
and the second operator corresponds to the target. Now 
consider the bit-wise application of seven CNOT gates 
between logically encoded blocks of data [Fig. 12 . First 
the stabilizer set must remain invariant, i.e, 



g = {K l ® K j } -> {IT ® K 3 } V 



(70) 



B. Two-qubit gate. 

A two-qubit logical CNOT operation can also be ap- 
plied in the same transversal way. For un-encoded qubits, 
a CNOT operation performs the following mapping on 



Table [V] details the transformation for all the stabilizers 
under seven bit-wise CNOT gates, demonstrating that 
this operation is valid on the [[7, 1,3]] code. The trans- 
formations in Eq. [69] are trivially extended to the logical 
space, showing that seven bit- wise CNOT gates invoke a 
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logical CNOT operation. 

i® z ^ z® z, 



Z®I- 
I®X 



Z®I, 
■I®X. 



(71) 



The issue of Fault-tolerance with these logical opera- 
tions should be clear. The X,Z, H and P gates are 
trivially fault-tolerant since the logical operation is per- 
formed through seven bit- wise single qubit gates. The 
logical CNOT is also fault-tolerant since each two-qubit 
gate only operates between counterpart qubits in each 
logical block. Hence if any gate is inaccurate, then at 
most a single error will be introduced in each block. 



7 



i 



FIG. 12 Bit-wise application of a CNOT gate between two 
logical qubits. Since each CNOT only couples correspond- 
ing qubits in each block, this operation is inherently fault- 
tolerant. 



In contrast to the [[7,1,3]] code, let us also take a quick 
look at the [[5,1,3]] code. As mentioned in section |VIII 
the [[5,1,3]] 



code is a non-CSS code, meaning the the 
Clifford group of gates cannot be fully implemented in a 
transversal manner. To see this clearly we can examine 
how the stabilizer group for the code transforms under a 
transversal Hadamard operation, 



/X Z Z X r 
I X Z Z X 
X I X Z Z 

K Z X I X z, 



'z x x z r 

I z X X z 

Z I z X X 

\X Z I Z Xj 



(72) 



The stabilizer group is not preserved under this transfor- 
mation, therefore the transversal Hadamard operation is 
not valid for the [[5,1,3]] code. One thing to briefly note is 
that there is a method for performing logical Hadamard 
and phase gates on the [[5,1,3]] code (Got97). However, it 
essentially involves performing a valid, transversal, three- 
qubit gate and then measuring out two of the logical an- 
cillae. 



While these gates are useful for operating on quantum 
data, they do not represent a universal set for quantum 
computation. In fact it has been shown that by using the 
stabilizer formalism, these operations can be efficiently 
simulated on a classical device (|Got98l IAG04p . In or- 
der to achieve universality one of the following gates are 
generally added to the available set, 



T = 



1 

e i7r / 4 



(73) 



or the Toffoli gate (|Tof8ip . However, neither of these 
two gates are members of the Clifford group and apply- 
ing them in a similar way to the other gates will trans- 
form the stabilizers out the group and consequently does 
not represent a valid operation. Circuits implementing 
these two gates in a fault-tolerant manner have been de- 
veloped (INC00I ICC991 ISl05l ISFH08L but at this stage 
the circuits are complicated and resource intensive. This 
has practical implications to encoded operations. If uni- 
versality is achieved by adding the T gate to the list, ar- 
bitrary single qubit rotations require long gate sequences 
(utilizing the Solovay-Kitaev theorem (|Kit97| IDN06P ) to 
approximate arbitrary logical qubit rotations and these 
sequences often require many T gates (Fow05). Finding 
more efficient methods to achieve universality on encoded 
data is therefore still an active area of research. 



XII. FAULT-TOLERANT CIRCUIT DESIGN FOR 
LOGICAL STATE PREPARATION 

Section [X] introduced the basic rules for fault-tolerant 
circuit design and how these rules lead to the thresh- 
old theorem for concatenated error correction. However, 
what does a full fault-tolerant quantum circuit look like? 
Here, we introduce a full fault-tolerant circuit to prepare 
the [[7, 1, 3]] logical |0) state. As the [[7, 1, 3]] code is a sin- 
gle error correcting code, we use the one-to-one definition 
of Fault-tolerance and therefore only need to consider the 
propagation of a single error during the preparation (any 
more that one error during correction represents a higher 
order effect and is ignored). 

As described in Section |VII| logical state preparation 
can be done by initializing an appropriate number of 
physical qubits and measuring each of the X stabiliz- 
ers that describe the code. Therefore, a circuit which 
allows the measurement of a Hermitian operator in a 
fault-tolerant manner needs to be constructed. The gen- 
eral structure of the circuit used was first developed by 
Shor (Sho96), however it should be noted that several 
more recent methods for fault-tolerant state preparation 
and correction now exist (|Ste97al ISte02l IDA07[) which 
are more efficient than Shor's original method. 

The circuits shown in Fig. |13| a, and |13b , which measure 
the stabilizer K 1 = IIIXXXX are not fault-tolerant, 
since a single ancilla is used to control each of the four 
CNOT gates. Instead, four ancilla qubits are used which 
are prepared in the state \A) = (|0000) + |llll))/\/2. 
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TABLE V Transformations of the [[7, 1, 3]] stabilizer set under the gate operation [/ =CNOT' 87 , where 5 -> JJ^f/. Note that 
the transformation does not take any stabilizer outside the group generated by 7f l ® TP (i, j) £ [1, .., 6], therefore (7 =CNOT® 7 
represents a valid operation on the codespace. 



This can be done by initializing four qubits in the |0) 
state and applying a Hadamard then a sequence of CNOT 
gates. Each of these four ancilla are used to control a 
separate CNOT gate, after which the ancilla state is de- 
coded and measured. By ensuring that each CNOT is 
controlled via a separate ancilla, any X error will only 
propagate to a single qubit in the data block. However, 
during the preparation of the ancilla state there is the 
possibility that a single X error can propagate to mul- 
tiple ancilla, which arc then fed forward into the data 
block. In order to combat this, the ancilla block needs to 
be verified against possible X errors. Tracking through 
all the possible locations where a single X error can occur 
during ancilla preparation leads to the following unique 
states. 
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(74) 



From these possibilities, the last four states have a dif- 
ferent parity between the first and forth qubit. Hence to 
verify this state, a fifth ancilla is added, initialized and 
used to perform a parity check on the ancilla block. This 
fifth ancilla is then measured. If the result is |0), the an- 
cilla block is clean and can be coupled to the data. If the 
ancilla result is |1), then either a single error has occured 
in the ancilla preparation or on this verification qubit. In 
either case, the entire ancilla block is re-initialized and 
the ancilla prepared again. This is continued until the 
verification qubit is measured to be |0) [Fig. 14 . The 



re-preparation of the ancilla block protects against X er- 
rors, which can propagate forward through the CNOT 
gates. Z errors on the other hand, propagate in the other 
direction. Any Z error which occurs in the ancilla block 
will propagate straight through to the final measurement. 
This results in the measurement not corresponding to the 
eigenstate the data is projected to and can possibly result 
in mis-correction once all stabilizers have been measured. 



To protect against this, each stabilizer is measured 2-3 
times and a majority vote of the measurement results 
taken. As any additional error represents a second or- 
der process, if the first or second measurement has been 
corrupted by an induced Z error, then the third mea- 
surement will only contain additional errors if a higher 
order error process has occurred. Therefore, we are free 
to ignore this possibility and assume that the third mea- 
surement is error free. The full circuit for [[7,1,3]] state 
preparation is shown in Fig. |15[ where each stabilizer is 
measured 2-3 times. The total circuit requires a mini- 
mum of 12 qubits (7-data qubits and a 5-qubit ancilla 
block). 

As you can see, the circuit constructions for full fault- 
tolerant state preparation (and error correction) are not 
simple circuits. However, they are easy to design in 
generic ways when employing stabilizer coding. 



XIII. LOSS PROTECTION 

So far we have focused the discussion on correction 
techniques which assume that error processes maintain 
the assumption of a qubit structure to the Hilbert space. 
As we noted in section |III.C[ the loss of physical qubits 
within the computer violates this assumption and in 
general requires additional correction machinery beyond 
what we have already discussed. 

For the sake of completeness, this section examines 
some correction techniques for qubit loss. Specifically, 
we detail one such scheme which was developed with sin- 
gle photon based architectures in mind. 

Protecting against qubit loss requires a different ap- 
proach than other general forms of quantum errors such 
as environmental decoherence or systematic control im- 
perfections. The cumbersome aspect related to correct- 
ing qubit loss is detecting the presence of a qubit at the 
physical level. The specific machinery that is required for 
loss detection is dependent on the underlying physical ar- 
chitecture, but the basic principal is that the presence or 
absence of the physical qubit must be determined with- 
out discriminating the actual quantum state. 

Certain systems allow for loss detection is a more 
convenient way than others. Electronic spin qubits, 
for example, can employ Single Electron Transistors 
(SET) to detect the presence or absence of the charge 
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FIG. 13 Three circuits which measure the stabilizer K 1 . Fig a) represents a generic operator measurement where a multi-qubit 
controlled gate is available. Fig. b) decomposes this into single- and two-qubit gates, but in a non-fault-tolerant manner. Fig. 
c) introduces four ancilla such that each CNOT is controlled via a separate qubit. This ensures Fault-tolerance. 



without performing measur ement on the spin degree of 
freedom (IDSOOl IOCJH05I lAJW+Oll) . Optics in con- 
trast requires more complicated non-demolition measure- 
ment (|MW84I IIHY85I lPOW+041 IMNBS05P - This is due 
to the fact that typical photonic measurement is per- 
formed via photo-detectors which have the disadvantage 
of physically destroying the photon. 

Once the detection of the presence of the physical qubit 
has been performed, a freshly initialized qubit can be 
injected to replace the lost qubit. Once this has been 
completed, the standard error correcting procedure can 
correct for the error. A freshly initialized qubit state, 
|0) can be represented as projective collapse of a general 
qubit state, \ip), as, 



ex \i>) + Z\iP) 



(75) 



If we consider this qubit as part of an encoded block, 
then the above corresponds to a 50% error probability of 
experiencing a phase flip on this qubit. Therefore, a loss 
event that is corrected by non-demolition detection and 



standard QEC essentially guarantees a correction event 
in the QEC cycle. Therefore the probability of loss needs 
to be at a comparable rate to standard errors as the cor- 
rection cycle after a loss detection event will, with high 
probability, detect and correct the error. 

Additionally, if a loss event is detected and the qubit 
replaced, the error detection code shown in section VI 
becomes a single qubit correction code. This is due to 
the fact that erasure errors have known locations. Conse- 
quently error detection is sufficient to perform full correc- 
tion, in contrast to non-erasure errors where the location 
is unknown. 

A second method for loss correction is related to sys- 
tems that have high loss rates compared to systematic 
and environmental errors. The most prevalent in op- 
tical systems. Due to the high mobility of single pho- 
tons and their relative immunity to environmental in- 
teractions, loss is a major error channel that generally 
dominates over other error sources. The use of error 
detection and correction codes for photon loss is unde- 
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FIG. 14 Circuit required to measure the stabilizer K , fault-tolerantly. A four qubit GHZ state is used as ancilla with the 
state requiring verification against multiple X errors. After the state has passed verification it is coupled to the data block and 
a syndrome is extracted. 



sirable due to the need for non-demolition detection of 
the lost qubit. While techniques exist for measuring the 
presence or absence of a photon without direct det ection 
have been developed and implemented (jPOW + 04[) . they 
require multiple ancilla photons and controlled interac- 
tions. Ultimately it is more desirable to redesign the loss 
correction code such that it can be employed directly 
with photo-detection rather than more complicated non- 
demolition techniques. 



One such scheme was developed by Ralph, Hayes and 
Gilchrist in 2005 (RHG05). This scheme was a more 
efficient extension of an original Parity encoding method 
developed by Knill, Laflamme and Milburn to protect 
against photon loss in their controlled-a^ gate i'KLMOl j). 
The general Parity encoding for a logical qubit is an N 
photon GHZ state in the conjugate basis, i.e, 
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where |±) = (|0) ± \1))/V2. The motivation with this 
type of encoding is that measuring any qubit in the |0, 1) 
basis simply removes it from the state, reducing the re- 



sulting state by one, i.e, 
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(77) 

where Po,i.jv are the projectors corresponding to mea- 
surement in the |0, 1) basis on the N th qubit (up to nor- 
malization). The effect for the \1) L state is similar. Mea- 
suring the N th qubit in the |0) state simply removes it 
from the encoded state, reducing the logical zero state 
by one, while measuring the N th qubit as |1) enacts a 
logical bit flip at the same time as reducing the size of 
the logical state. However, since the measurement result 
is known, this encoded bit flip can be corrected for. 

Instead of introducing the full scheme developed 
in (|RHG05p . we instead just give the general idea of 
how such encoding allows for loss detection without non- 
demolition measurements. Photon loss in this model 
is assumed equivalent to measuring the photon in the 
|0),|1) basis, but not knowing the answer [Sec III.C 



Our ignorance of the measurement result could lead to a 
logical bit flip error on the encoded state, therefore we 
require the ability to protect against logical bit flip er- 
rors on the above states. As already shown, the 3-qubit 
code allows us to achieve such correction. Therefore the 
final step in this scheme is encoding the above states into 
a redundancy code (a generalized version of the 3-qubit 
code), where an arbitrary logical state, \ip) l is now given 
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FIG. 15 Circuit required to prepare the [[7, 1,3]] logical |0) state fault-tolerantly. Each of the X stabilizers are sequentially 
measured using the circuit in Fig. |14| To maintain Fault-tolerance, each stabilizer is measured 2-3 times with a majority vote 
taken. 
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A. Subsystem Codes 



= a|0)f|0)^...|0)f + /3|l)f|l)f. 



(78) 



where 10)^ ,\l) N are the parity encoded states shown in 
Eq. [77] and the fully encoded state is ^-blocks of these 
parity states. 

This form of encoding protects against the loss of 
qubits by first encoding the system into a code structure 
that allows for the removal of qubits without destroy- 
ing the computational state and then protecting against 
logical errors that are induced by loss events. In effect 
it maps errors that are un-correctable by standard QEC 
to error channels that are correctable, in this case qubit 
loss — > qubit bit-flips. This is common with pathological 
error channels. If a specific type of error violates the stan- 
dard "qubit" assumption of QEC, additional correction 
techniques are always required to map this type of error 
to a correctable form, consequently additional physical 
resources are usually needed. 



XIV. SOME MODERN DEVELOPMENTS IN ERROR 
CORRECTION 

Up until this stage we have restricted our discussions 
on error correction to the most basic principals and codes. 
The ideas and methodologies we have detailed represent 
the introductory techniques that were developed when 
error correction was first proposed. For those readers 
who are only looking for a basic introduction to the field, 
you can quite easily skip the remainder of this paper. 

Providing a fair and encompassing review of the more 
modern and advanced error correction techniques that 
have been developed is far outside our goal for this re- 
view. However, we would be remiss if we did not briefly 
examine some of the more advanced error correction tech- 
niques that have been proposed for large scale quantum 
information processing. For the remainder of this discus- 
sion we choose two closely related error correction tech- 
niques, subsystem coding and topological coding which 
have been receiving significant attention in the fields of 
architecture design and large scale quantum information 
processing. While some readers may disagree, we review 
these two modern error correction protocols because they 
are currently two of the most useful correction techniques 
when discussing the physical construction of a quantum 
computer. 

We again attempt to keep the discussion of these tech- 
niques light and provide specific examples when possible. 
However, it should be stressed that these error correct- 
ing protocols are far more complicated than the basic 
codes shown earlier. Topological error correction alone 
has since its introduction, essentially become its own 
research topic within the broader error correction field. 
Hence we encourage the reader who is interested to refer 
to the cited articles below for more rigorous and detailed 
treatment of these techniques. 



Quantum subsystem codes (!Bac06p are one of the 
newer and highly flexible techniques to implement quan- 
tum error correction. The traditional stabilizer codes 
that we have reviewed are more formally identified as 
subspace codes, where information is encoded in a rele- 
vant coding subspace of a larger multi-qubit system. In 
contrast, subsystem coding identifies multiple subspaces 
of the multi-qubit system as equivalent for storing quan- 
tum information. Specifically, multiple states are identi- 
fied with the logical \0) L and |1) L states. 

The primary benefit to utilizing subsystem codes is 
the general nature of their construction. Moving from 
smaller to larger error correcting codes is conceptually 
straightforward, error correction circuits are much sim- 
pler to construct when encoding information in multiple 
subsystems (|AC07) and the generality of their construc- 
tion introduces the ability to perform dynamical code 
switching in a fault-tolerant manner (SEDH07). This fi- 
nal property gives subsystem coding significant flexibility 
as the strength of error correction within a quantum com- 
puter can be changed, fault-tolerantly, during operation 
of the device. 

As with the other codes presented in this review, sub- 
system codes are stabilizer codes but now defined over 
a square lattice. The lattice dimensions represent the 
X and Z error correction properties and the size of the 
lattice in cither of these two dimensions dictates the to- 
tal number of errors the code can correct. In general, a 
C (77-1,712) subsystem code is defined over a n\ x n 2 square 
lattice which encodes one logical qubit into n \ni physical 
qubits with the ability to correct at least j Z errors 

and at least X errors. Again, keeping with the 

spirit of this review, we instead focus on a specific exam- 
ple, the C(3,3) subsystem code. This code, encoding 9 
physical qubits into one logical qubit can correct for one 
X and one Z error. In order to define the code structure 
we begin with a 3 x 3 lattice of qubits, where each qubit 
is identified with the vertices of the lattice (note that 
this 2D structure represents the structure of the code, it 
does not imply that a physical array of qubits must be 
arranged into a 2D lattice) . Fig. [16] illustrate three sets 
of stabilizer operators which are defined over the lattice. 
The first group, illustrated in Fig. [I6J1. is the stabilizer 
group, S, which is generated by the operators, 



s = (X^X. 



i+ly 



ieZ 2 ;jeZ 2 ), (79) 



where we retain the notation utilized in (|AC071 SEI 
Ui^ and U*.j represent an operator, U, acting on all 
qubits in a given row, i, or column, j, respectively, and 
Z2 = {1,2}. The second relevant subsystem is known as 
the gauge group [Fig. 16 d.], T, and is described via the 



non-Abelian group generated by the pairwise operators 



T ={XijX i+ i tj I i eZ 2 ;j e Z 3 ), 
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FIG. 16 Stabilizer structure for the C(3,3) code. Fig a. gives 
two of the four stabilizers from the group S. Fig. b. illustrates 
one of the four encoded Pauli operators from each subsystem 
defined with the Gauge group, T. Fig. c. gives the two logical 
operators from the group C which enact valid operations on 
the encoded qubit. 



The third relevant subsystem is the logical 
space [Fig. 16;], £, which can be defined through 
the logical Pauli operators 



l,*/ 



(81) 



which when combined with S form a non-Abelian group. 

The stabilizer group S, defines all relevant code states, 
i.e. every valid logical space is a +1 eigenvalue of this 
set. For the C(3,3) code, there are a total of nine physi- 
cal qubits and a total of four independent stabilizers in S, 
hence there are five degrees of freedom left in the system 
which can house 2 5 logical states which are simultaneous 
eigenstates of S. This is where the gauge group, T, be- 
comes relevant. As the gauge group is non-Abelian, there 
is no valid code state which is a simultaneous eigenstate 
of all operators in T. However, if you examine closely 
there are a total of four encoded Pauli operations within 
T. 



Fig 16 



illustrates two such operators. As all ele- 
ments of T commute with all elements of S we can iden- 
tify each of these four sets of valid "logical" qubits to 
be equivalent, i.e. we define {|0) L , |1) L } pairs which are 
eigenstates of S and an abelian subgroup of T and then 
ignore exactly what gauge group we are in (each of the 
four possible |0) L states can be used to store a single log- 



ical qubit in the |0) state, regardless of which particular 
|0) L gauge state we are in). Hence, each of these gauge 
states represent a subsystem of the code, with each sub- 
system logically equivalent. 

The final group we considered is the logical group C. 
This is the set of two Pauli operators which enact a logical 
X or Z gate on the encoded qubit regardless of the gauge 
choice and consequently represent true logical operations 
to our encoded space. 

In a more formal sense, the definition of these three 
group structures allows us to decompose the Hilbert 
space of the system. If we let H denote the Hilbert space 
of the physical system, S forms an Abelian group and 
hence can act as a stabilizer set denoting subspaces of H. 
If we describe each of these subspaces by the binary vec- 
tor, e, formed from the eigenvalues of the stabilizers, S, 
then each subspace splits into a tensor product structure 
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where elements of T act only on the subsystem Tlj- and 
the operators C act only on the subsystem He- There- 
fore, in the context of storing qubit information, a logical 
qubit is encoded into the two dimensional subsystem Tic ■ 
As the system is already stabilized by operators in S and 
the operators in T act only on the space Hr, qubit in- 
formation is only altered when operators in the group C 
act on the system. 

This formal definition of how subsystem coding works 
may be more complicated than the standard stabilizer 
codes shown earlier, but this slightly more complicated 
coding structure has significant benefits when we consider 
how error correction is performed. In general, to perform 
error correction, each of the stabilizers of the codespace 
must be checked to determine which eigenvalue changes 
have occurred due to errors. In the case of subsystem 
code this would appear to be problematic. The stabilizer 
group, S ', consist of qubit operators that scale with the 
size of the code. In our specific example, each of the X 
and Z stabilizers are six-dimensional (and in general, for 
a ni x n 2 lattice, the X stabilizers are 2ri\ dimensional 
and the Z stabilizers are 2ri2 dimensional) . If techniques 
such as Shor's method [Section XII were used, we would 



need to prepare a large ancilla state to perform fault- 
tolerant correction, which would also scale linearly with 
the size of the code, this is clearly undesirable. However, 
due to the gauge structure of subsystem codes we are able 
to decompose the error correction procedure ()AC07p . 

Each of the stabilizers in S are simply the product of 
certain elements from T, for example, 



Xi, 1X1,2X1, 3X2, 1X2, 2X2, 3 € S 

= (^1, 1^2, l)- ^1,2X2, 2) -(Xi, 3X2, 3) 



ST. 



(83) 



Therefore if we check the eigenvalues of the three, 2- 
dimensional operators from T we are able to calculate 
what the eigenvalue is for the 6-dimensional stabilizer. 
This decomposition of the stabilizer set for the code can 
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only occur since the decomposition is in terms of opera- 
tors from T which, when measured, has no effect on the 
logical information encoded within the system. In fact, 
when error correction is performed the gauge state of the 
system will almost always change based on the order in 
which the eigenvalues of the gauge operators are checked. 

This exploitation of the gauge properties of subsystem 
coding is extremely beneficial for fault-tolerant designs 
for correction circuits. As the stabilizer operators can 
now be decomposed into multiple 2-dimensional opera- 
tors, fault-tolerant circuits for error correction do not 
require any encoded ancilla states. Furthermore, if we 
decide to scale the code-space to correct more errors (in- 
creasing the lattice size representing the code) we do 
not require measuring operators with higher dimension- 
ality. Fig 17 taken from Ref. (|AC07[1 illustrates the fault- 



tolerant circuit constructions for Bacon-Shor subsystem 
codes. As each ancilla qubit is only coupled to two data 



(a) 



(b) 



j+i, k 



1+) 



k,j+l 



H 



trates circuits to perform fault-tolerant switching from 
the C(3,3) and C(3,5) subsystem code. As noted before, 
the C(3,3) is a single X, single Z error correcting code 
while the C(3,5) is a single X, two Z error correcting 
code. We will not detail why these circuits successfully 
implement fault-tolerant code switching, instead we en- 
courage readers to refer to Ref. ([SEDH07P for further 
details. 



2,j 



3,j 



|0> 

|0> 



H 



\0)~H 



2,3 



3,j 



i,3 



5,3 



FIG. 18 (From Ref. (SEDH07)). Circuit to convert from the 
C(3,3) subsystem code to the C(5,3) code for one column, j, 
of the lattice structure of the code. 



FIG. 17 (From Ref. ljAC07J0 Circuits for measuring the 
gauge operators and hence performing error correction for 
subsystem codes. Fig. a. measures, fault-tolerantly, the op- 
erator Xj^Xj + \^ with only one ancilla. Fig. b. measures 
Zk,jZk,j+i- The results of these two qubit parity checks can 
be used to calculate the parity of the higher dimensional sta- 
bilizer operators of the code. 

qubits, no further circuit constructions are required to 
ensure fault-tolerance. The classical results from these 
2-dimensional parity checks are then combined to calcu- 
late the parity of the higher dimensional stabilizer of the 
subsystem code. 

A second benefit to utilizing subsystem codes is the 
ability to construct fault-tolerant circuits to perform dy- 
namical code switching. When using more traditional 
error correction codes it is difficult, if not impossible, to 
fault-tolerantly switch between codes with different error 
correcting properties. The Steane [[7,1,3]] code is a single 
error correcting code for both X and Z channels. If dur- 
ing the operation of a quantum computer, the user wished 
to increase the error correcting power of their code to two 
errors in the X and Z channel they would first decode 
the quantum data and then re-encode with the higher 
distance code. This is clearly a non fault-tolerant proce- 
dure as any error occurring on the decoded information 
will cause catastrophic failure. Due to the general lat- 
tice structure of subsystem codes switching too and from 
higher order codes can be achieved without decoding and 
re-encoding information, allowing the user of the com- 
puter to dynamically adjust the error correction during 
the computation. 



Figs. Ft8] [19] and |20j taken from Ref. (|SEDH07|) illus- 
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FIG. 19 (From Ref. I|SEDH07|) '). Downconversion from the 
C(5,3) code to the C(3,3)code . V is the gate sequence in 

Fig.^rn 
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FIG. 20 (From Ref. (SEDH07)) X parity measurement under 
C(5,3) for one row, i, of the lattice structure. 
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B. Topological Codes 

A similar coding technique to the Bacon-Shor subsys- 
tem codes is the idea of topological error correction, first 
introduced with the Toric code of Kitaev in 1997 (|Kit97[l . 
Topological coding is similar to subsystem codes in that 
the code structure is defined on a lattice (which, in gen- 
eral, can be be of dimension > 2) and the scaling of the 
code to correct more errors is conceptually straightfor- 
ward. However, in topological coding schemes the pro- 
tection afforded to logical information relies on the un- 
likely application of error chains which define non-trivial 
topological paths over the code surface. 

Topological error correction is a complicated area of 
quantum error correction and fault-tolerance and any at- 
tempt to fairly summarize the field is not possible within 
this review. In brief, there are two ways of approach- 
ing the problem. The first is simply to treat topologi- 
cal codes as a class of stabilizer codes over a qubit sys- 
tem. This approach is more amenable to current cur- 
rent information technologies and is being adapted to 
metho ds in clus ter state computing (RHGOT] IFG08|I . 
optics (|DFS+08l IDMN08I), io n-traps (|SJ08)) and super- 
conducting systems ( IFI + 02[| . The second approach is 
to construct a physical Hamiltonian model based on 
the structure of the topological code. This leads to 
the more complicated field on anyonic quantum com- 
putation (Kit97). By translating a coding structure 
into a physical Hamiltonian system, excitations from the 
ground state of this Hamiltonian exhibit natural robust- 
ness against local errors (since their physical Hamilto- 
nian symmetries reflect the coding structure imposed). 
Specifically, quasi-particles arising from a Hamiltonian 
approach to quantum codes exhibit fractional quantum 
statistics (they acquire fractional phase shifts when their 
positions are exchanged twice with other anyons, in 
contrast to Bosons or Fermions which always acquire 
±1 phase shifts). The unique properties of anyons 
therefore allow for natural, robust, error protection and 
anyon/anyon interactions are performed by rotating any- 
ones around each other. However, the major issue with 
this model is that it relies on quasi-particle excitations 
that do not, in general, arise naturally. Although certain 
physical systems have been shown to exhibit anyonic ex- 
citations, most notably in the fractional quantum hall ef- 
fect (|NSS + 08|l the ability to first manufacture a reliable 
anyonic system in addition to reliably design and con- 
struct a large scale computing system based on anyons is 
a daunting task. 

As the re are sev eral extremely good discussions of both 
anyonic (|NSS+08|l and non-anyonic topological comput- 
ing (|DKLP02l IFSG08I IFG08|) we will not review any 
of the anyonic methods for topological computing and 
simply provide a brief example of one topological cod- 
ing scheme, namely the surface code (BK01; DKLP02 
IFSG08|1 . The surface code for quantum error correction 
is an extremely good error correction model for several 
reasons. As it is defined over a 2-dimensional lattice of 



qubits it can be implemented on architectures that only 
allow for the coupling of nearest neighbor qubits (rather 
than the arbitrary long distance coupling of qubits in 
separate regions of the computer). The surface code also 
exhibits one of the highest fault-tolerant thresholds of 
any quantum error correction scheme, recent simulations 
estimate a threshold approaching 1% (RHG07). Finally, 
the surface code can correct problematic error channels 
such as qubit loss and leakage naturally. 

The surface code, as with subsystem codes, is a sta- 
bilizer code defined over a 2-dimensional qubit lattice. 
Fig. [2T| illustrates. We now identify each edge of the 2D 
lattice with a physical qubit. The stabilizer set consists 
of two types of operators, the first is the set of iT® 4 oper- 
ators which circle every lattice face (or plaquette). The 
second is the set of X® 4 operators which encircle every 
vertex of the lattice. The stabilizer set is consequently 
generated by the operators, 



A p — 



B„ — 



X, 



(84) 



jeb(p) 



where b(p) is the four qubits surrounding a plaquette and 
s(v) is the four qubits surrounding each vertex in the 
lattice and identity operators on the other qubits are im- 
plied. First note that all of these operators commute 
as any plaquette and vertex stabilizer will share either 
zero or two qubits. If the lattice is not periodic in ei- 
ther dimension, this stabilizer set completely specifies one 
unique state, i.e. for a N x N lattice there are 27V 2 qubits 
and 27V 2 stabilizer generators. Hence this stabilizer set 
defines a unique multi-qubit entangled state which is gen- 
erally referred to as a ■■clean" surface. Detailing exactly 
how this surface can be utilized to perform robust quan- 
tum computation is far outside the scope of this review 
and there are several papers to which such a discussion 
can be referred (IRH07I IPJTG071 IFSG08I IFG081) . Instead, 
we can quite adequately show how robust error correction 
is possible by simply examining how a "clean" surface can 
be maintained in the presence of errors. The X and Z 
stabilizer sets, A p and B v define two equivalent 2D lat- 
tices which are interlaced, Fig. 22 illustrates. If the total 
2D lattice is shifted along the diagonal by half a cell then 
the operators B v are now arranged around a plaquette 
and the operators A p are arranged around a lattice ver- 
tex. Since protection against X errors are achieved by 
detecting eigenvalue flips of Z stabilizers and visa-versa, 
these two interlaced lattices correspond to error correc- 
tion against X and Z errors respectively. Therefore we 
can quite happily restrict our discussion to one possible 
error channel, for example correcting X errors (since the 
correction for Z errors proceeds identically when consid- 
ering the stabilizers B v instead of A p ) . 

Fig [23] i. illustrates the effect that a singe X error 
has on a pair of adjacent plaquettes. Since X and Z 
anti-commute, a single bit-flip error on one qubit in the 
surface will flip the eigenvalue of the Z® 4 stabilizers on 
the two plaquettes adjacent to the respective qubit. As 
single qubit errors act to flip the eigenvalue of adjacent 
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Qubits 




FIG. 21 General structure of the surface code. The edges 
of the lattice correspond to physical qubits. The four qubits 
surrounding each face (or plaquette) are +1 eigenstates of the 
operators A p while the four qubits surrounding each vertex 
are +1 eigenstates of the operators B v . If all eigenstate con- 
ditions are met, a unique multi-qubit state is defined as a 
"clean" surface. 



plaquette stabilizers we examine how chains of errors af- 
fect the surface. Figs 23 3. and Fig. 23 examine two 



longer chains of errors. As you can sec, if multiple errors 
occur, only the eigenvalues of the stabilizers associated 
with the ends of the error chains flip. Each plaquette 
along the chain will always have two X errors occurring 
on different boundaries and consequently the eigenvalue 
of the Z® A stabilizer around these plaquettes will flip 
twice. 

If we now consider an additional ancilla qubit which 
sits in the center of each plaquette and can couple to the 
four surrounding qubits, we can check the parity by run- 
ning the simple parity circuit shown in Fig |24j If we as- 
sume that we initially prepare a perfect "clean" surface, 
we then, at some later time, check the parity of every 
plaquette over the surface. If X errors have occurred on 
a certain subset of qubits, the parity associated with the 
endpoints of error chains will have flipped. We now take 
this 2-dimensional classical data tree of eigenvalue flips 
and pair them up into the most likely set of error chains. 
Since it is assumed that the probability of error on any 
individual qubit is low, the most likely set of errors which 
reflects the eigenvalue changes observed is the minimum 
weight set (i.e. connect up all plaquettes where eigenval- 
ues have changed into pairs such that the total length of 
all connections is minimized) . This classical data process- 
ing is quite common in computer science and minimum 
weight matching algorithms such as the Blossom pack- 
age (CR99 Kol08) have a running time polynomial in 



the total number of data points in the classical set. Once 
this minimal matching is achieved, we can identify the 
likely error chains corresponding to the end points and 
correction can be applied accordingly. 

The failure of this code is therefore dictated by error 
chains that cannot be detected through changes in pla- 
quette eigenvalues. If you examine Fig |25[ we consider 
an error chain that connects one edge of the surface lat- 
tice to another. In this case every plaquette has two 
associated qubits that have experienced a bit flip and no 
eigenvalues in the surface have changed. Since we have 
assumed that we are only wishing to maintain a "clean" 
surface, these error chains have no effect, but when one 
considers the case of storing information in the lattice, 
these types of error chains correspond to logical errors 
on the qubit (FSG08). Hence undetectable errors are 
chains which connect boundaries of the surface to other 
boundaries (in the case of information processing, qubits 
are artificial boundaries within the larger lattice surface) 
It should be stressed that this is a simplified description 
of the full protocol, but it does encapsulate the basic idea. 
The important thing to realize is that the failure rate of 
the error correction procedure is suppressed, exponen- 
tially with the size of the lattice. In order for a series of 
single qubit errors to be undetectable, they must form a 
chain connecting one boundary in the surface with an- 
other. If we consider an error model where each qubit 
experiences a bit flip, independently, with probability p, 
then an error chain of one occurs with probability p, error 
chains of weight two occur with probability 0(p 2 ), chains 
of three 0(p 3 ) etc... If we have an N x N lattice and we 
extend the surface by one plaquette in each dimension, 
then the probability of having an error chain connecting 
two boundaries will drop by a factor of p 2 (two extra 
qubits have to experience an error one on each bound- 
ary). Extending and N x N lattice by one plaquette in 
each dimension requires O(N) extra qubits, hence this 
type of error correcting code suppresses the probability 
of having undetectable errors exponentially with a qubit 
resource cost which grows linearly. 

As we showed in Section [X] standard concatenated 
coding techniques allow for a error rate suppression which 
scales with the concatenation level as a double exponen- 
tial while the resource increase scales exponentially. For 
the surface code, the error rate suppression scales ex- 
ponentially while the resource increase scales linearly. 
While these scaling relations might be mathematically 
equivalent, the surface code offers much more flexibility 
at the architectural level. Being able to increase the error 
protection in the computer with only a linear change in 
the number of physical qubits is far more beneficial than 
using an exponential increase in resources when utilizing 
concatenated correction. Specifically, consider the case 
where a error protected computer is operating at a logi- 
cal error rate which is just above what is required for an 
algorithm. If concatenated error correction is employed, 
then adding another later of correction will not only in- 
crease the number of qubits by an exponential amount, 
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FIG. 22 The surface code imbeds two self similar lattices that are interlaced, generally referred to as the primal and dual 
lattice. Fig. a. illustrates one lattice where plaquettes are defined with the stabilizers A p . Fig b. illustrates the dual structure 
where plaquettes are now defined by the stabilizer set B v . The two lattice structures are interlaced and are related by shifting 
along the diagonal by half a lattice cell. Each of these equivalent lattices are independently responsible for X and Z error 
correction 
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FIG. 23 Examples of error chains and their effect on the eigenvalues for each plaquette stabilizer, a). A single X error causes 
the parity of two adjacent cells to flip, b) and c). Longer chains of errors only cause the end cells to flip eigenvalue as each 
intermediate cell will have two X errors and hence the eigenvalue for the stabilizer will flip twice. 



but it will also drop the effective logical error rate far 
below what is actually required. In contrast, if surface 
codes are employed, we increase the qubit resources by 
a linear factor and drop the logical error rate sufficiently 
for successful application of the algorithm. 



We now leave the discussion regarding topological cor- 
rection models. We emphasize again that this was a 
very broad overview of the general concept of topolog- 
ical codes. There are many details and subtleties that 
we have deliberately left out if this discussion and we 
urge the reader, if they are interested, to refer to the ref- 
erenced articles for a much more thorough treatment of 
this topic. 



XV. CONCLUSIONS AND FUTURE OUTLOOK 

This review has hopefully provided a basic introduc- 
tion to some of the most important theoretical aspects 
of quantum error correction and fault-tolerant quantum 
computation. The ultimate goal of this discussion was 
not to provide a rigorous theoretical framework for QEC 
and fault-tolerance, but instead attempted to illustrate 
most of the important rules, results and techniques that 
have evolved out of this field. 

We not only covered the basic aspects of QEC through 
specific examples, but also we have briefly discussed how 
physical errors influence quantum computation and how 
these processes are interpreted within the context of 
QEC. One of the more important aspects of this review is 
the discussion related to the stabilizer formalism, circuit 
synthesis and fault-tolerant circuit construction. Stabi- 
lizers are arguably the most useful theoretical formalism 
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FIG. 24 a). Lattice structure to check the parity of a surface 
plaquette. An additional ancilla qubit is coupled to the four 
neighboring qubits that comprise each plaquette. b). Quan- 
tum circuit to check the parity of the Z® 4 stabilizer for each 
surface plaquette. 




FIG. 25 Example of a chain of errors which do not cause any 
eigenvalue changes in the surface. If errors connect bound- 
aries to other boundaries, the error correction protocol will 
not detect them. In the case of a "clean" surface, these error 
chains are invariants of the surface code. When computa- 
tion is considered, qubit information are artificial boundaries 
within the surface. Hence if error chains connect these infor- 
mation qubits to other boundaries, logical errors occur. 



in QEC as once it is sufficiently understood, most of the 
important properties of error correcting codes can be in- 
vestigated and understood largely by inspection. 

The study of quantum error correction and fault- 
tolerance is still and active area of QIP research. Al- 
though the library of quantum codes and error correction 
techniques are vast there is still a significant disconnect 
between the abstract framework of quantum coding and 
the more physically realistic implementation of error cor- 



rection for large scale quantum information. 

There are several future possibilities for the direction of 
quantum information processing. Even with the develop- 
ment of many of these advanced techniques, the physical 
construction and accuracy of current qubit fabrication is 
still insufficient to obtain any benefit from QEC. Many 
in the field now acknowledge that the future development 
of quantum computation will most likely split into two 
broad categories. The first is arguably the more physi- 
cally realistic, namely few qubit application in quantum 
simulation. 

Quantum simulation, i.e. using quantum systems to 
efficiently simulate other quantum systems was proposed 
by Richard Feynmann in the early 1980's (Fey82j and 
was one of the primary motivations to the development 
of the field. In the ideal case, it is argued that having 
access to a quantum computer with on the order of 10- 
100 physical qubits could allow for simulating physical 
systems large enough to be impractical for current classi- 
cal computers. If we limit our quantum array to the 100 
qubit level, then even implementing active error correc- 
tion techniques would not be desirable. Instead, higher 
quality fabrication and control as well as techniques in 
error avoidance (which require far less resources than er- 
ror correction) would instead be used in order to lower 
effective error rates below what is required to run few 
qubit applications. 

Beyond few qubit quantum simulation we move to tru- 
ely large scale quantum computation, i.e. implementing 
large algorithms such as Shor on qubit arrays well beyond 
1000 physical qubits. This would undoubtably require 
active techniques in error correction. Future work needs 
to focus on adapting the many codes and fault-tolerant 
techniques to the architectural level. As we noted in sec- 
tion |X.B[ the implementation of QEC at the design level 
largely influences the fault-tolerant threshold exhibited 
by the code itself. Being able to efficiently incorporate 
both the actual quantum code and the error correction 
procedures at the physical level is extremely important 
when developing an experimentally viable, large scale 
quantum computer. 

There are many differing opinions within the quan- 
tum computing community as to the future prospects 
for quantum information processing. Many remain pes- 
simistic regarding the development of a million qubit de- 
vice and instead look towards quantum simulation in the 
absence of active error correction as the realistic goal of 
quantum information. However, in the past few years, 
the theoretical advances in error correction and the fan- 
tastic speed in the experimental development of few qubit 
devices continues to offer hope for the near-term con- 
struction of a large scale device, incorporating many of 
the ideas presented within this review. While we could 
never foresee the possible successes or failures in quan- 
tum information science, we remain hopeful that a large 
scale quantum computer is still a goal worth pursuing. 
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